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PRECEDING rhSttL NOT FILMED. 


SUMMARY 


On June 18, 1976, a Scout D rocket was launched from Wallops Island, Virginia, 
carrying an atomic hydrogen maser oscillator system as the payload. The frequency 
of signals from the oscillator was monitored on the ground at Merritt Island, Florida, 
by using two hydrogen masers as comparison oscillators. The first-order doppler 
shift in the signals was eliminated by a go-return transponder link to the payload, and 
the resulting data, representing the relativistic shifts, were recovered and recorded. 
The objective of this experiment was to measure directly the effect of the gravitational 
potential on the frequency of an atomic hydrogen maser assuming it to be a ''proper” 
clock. 

In the experiment, a gravitational effect amounting to some 4. 5 parts in 10^^ was 

14 

measured with an oscillator having a stability better than 1 part in 10 . Therefore, 

to make the best possible use of the oscillator, we must account for all frequency shifts 
-15 

at the 2 to 5 X 10 level in Af/f in the system, and this, of course, includes all the 
phase variations that can cause such a shift to appear. 

This report presents the experiment, a description of the data now available and 
the manner in which they were processed, and the results. 

We wish to emphasize that this experiment was conducted by a large number of 
people and involved many different organizations. The program was managed by the 
George C. Marshall Space Flight Center of the National Aeronautics and Space Admin- 
istration (NASA), which was also responsible for designing and assembling the payload 
structure and testing it before launch. NASA's Goddard Space Flight Center performed 
the task of tracking and data acquisition, and NASA's Langley Research Center and the 
Vought Corporation were responsible for the vehicle. We gratefully acknowledge the 
assistant of the Jet Propulsion Laboratory in the solution of the probe's trajectory. 

The Smithsonian Astrc^hysical Observatory provided the scientific and teclinical leader- 
ship and constructed the probe and ground clocks and the doppler cancellation system.. 
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GRAVITATIONAL REDSHIFT SPACE-PROBE EXPERIMENT 


GP-A Final Report 
Contract NAS8-27969 


1. INTRODUCTION 

The test of the gravitational redshift predicted by Albert Einstein in 1907 had its 
beginnings with the development of atomic clocks in the 1950s; at that time, J. R. 
Zacharias' group at the Massachusetts Institute of Teclmology began working on 
cesium-beam clocks with an eye toward conducting an ejq)eriment between a clock 
placed on a mountain top and one in a valley to detect the predicted gravitational 
redshift. In a sense, Zacharias is a godfather of our present experiment. Its other 
godfather is N.F. Ramsey (1956), who discusses redshift measurement in Molecular 
Beams ; more importantly, however, it was Ramsey who invented the hydrogen maser, 
in 1959. 

The redshift test itself was first performed by R. V. Pound (Pound and Rebka, 

1960) using a 75-ft elevator shaft. His now classic \% test of the gravitational redshift 
was performed in I960 by using the then recently discovered Mossbauer effect. It is 
interesting to note that this terrestrial test was made over a height approximately the 
same as that of the Scout rocket used in the present e;q>eriment and discussed below. 
The test by Poimd and his coworkers, S.A. Rebka, Jr., and J. L. Snider, temporarily 
wiped out any further zeal to conduct a redshift experiment with clocks. To improve 
on their results would require far better clocks and enormously greater gravitational - 
potential differences between them. 

However, the space program, in fact, did evolve methods to obtain larger gravi- 
tational-potential differences, and the atomic hydrogen maser was invented in 1959. 

This new type of oscillator was developed into a stable clock that could resolve time 
-14 

differences at the 10 level. With the entire effect of the earth's gravity (consisting 
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of 7 parts in 10^®) made available, thanks to the space program of the National Aero- 
nautics and Space Administration (NASA), the prospect of a significant improvement 
over the 1% tests made by Pound became evident, and the redshift clock experiment 
was revived in 1966 by Ramsey, D. Kleppner, and R.F. C. Vessot (Ramsey et al . . 

1970). 

Unfortunately, circumstances favoring this experiment soon changed. The 
Vietnamese war turmoil at universities, the guns and butter disputes, and above all, 
the tremendously ejqienaive prospect of a Titan 3C vehicle and a 2000-lb payload all 
led to the termination of this version of the experiment. At about the same time, the 
maser group working at Varian Associates and Hewlett-Packard was successfully 
transplanted to the Smithsonian Astrophysical Observatory (SAO) to continue clock . 
development in support of satellite tracking, astronomy, and very long-baseline inter- 
ferometry. 

Redshift e^eriments with clocks remained dormant until 1970, when, at the Con- 
ference on Experimental Tests of Gravitation Theories at the California Institute of 
Technology, Vessot and M.W. Levine gave a paper (1971) describing a 24-hour eccentric- 
orbit redshift experiment and its capabilities. This talk prompted N.G. Roman and J. 
Mitchell to suggest an exceedingly eccentric orbit, essentially vertically up and down, 
achievable by the Scout rocket system. While this was a far more modest version than 
the original experiment, we were encouraged to pursue it. After some initial calculations, 
we concluded that it was indeed a viable possibility; the most serious limitation would be 
to keep the payload weight low enough for the Scout rocket to achieve both a reasonably 
high altitude and a sufficient time aloft. This necessitated a complete replanning of the 
(^erational strategy of the experiment. The payload system would have to recover 
very rapidly from the severe transitions in environment and the trauma of the laundi 
itself. The hazards of such a one-shot approach led to a very severe set of design 
constraints. 

This report describes the new departure by NASA into the field of relativity in the 
flying of a dedicated payload as a laboratory experiment to study gravitation. The 
laboratory itself consists of the earth and its gravity field and extends over the total 
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range of distance traveled by the probe vehicle and the earth station used to obtain the 
data, the two being connected by electromagnetic signals propagated to and from the 
probe. 
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2. THE EQUIVALENCE PRINCIPLE 


In 1907, Albert Einstein announced the principle of equivalence, which was doubt- 
less inspired by his reflections on the universality of free fall ** that all objects, regard 
less of their weight or composition, follow the same path when allowed to fall in a 
gravity field. Einstein extended his observation to state that "all freely falling, non- 
rotating laboratories are completely equivalent" and that all physics experiments done 
within the local confines of such laboratoiies are governed by the laws of special 
relativity. In particular, this means that it is impossible to distinguish, over small 
distances, the difference between a laboratory acted on by gravity and the same labora- 
tory accelerated by mechanical means. 

The equivalence principle is, in fact, the cornerstone of Einstein's General Theory 
of Relativity and of all other relativity theories that are based on a geometrical concept 
of space-time, the so-called metric theories. 

Metric theories of gravity - in particular. General Relativity, the most promising 
of all theories - arose from a merger of curved geometry with the four-dimensional 
space that Einstein and others used to describe physics in inertial (or unaccelerated) 
frames of reference. The description for such inertial frames Involves three dimen- 
sions of space plus the additional dimension of time and is called the Special Theory 
of Relativity. It is considered a flat space, since the four coordinate axes are straight 
lines and mathematically orthogonal. To include the effect of accelerations and gravi- 
tation, Einstein envisioned a warping of this coordinate system and used the geometry 
of curvilinear coordinates developed by Riemann and others to express the curvature. 
With extraordinary insight, he then related this curvature to physical parameters. 

The outcome is that physical experiments described in such a curvilinear coordinate 
frame include ^ aspects of acceleration and gravitation. Over small enough regions, 
the curvature is negligible and the "local" coordinates have the appearance of the 
flat space-time of special relativity. Figure 1 attempts to provide a visualization of 
this local connection between gravitation and acceleration. 
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"ALL LOCAL. FREELY FALLING, NONROTATING LABORATORIES ARE 
FULLY EQUIVALENT F(» THE PERFORMANCE OF ALL 

PHYSICAL EXPERIMENTS" A. EINSTEIN - 1907 - 


IT FOLLOWS THAT THERE IS MO WAY LOCALLY OF DISTINGUISHING 
BETWEEN THE PULL OF GRAVITY AND AN OPPOSITELY DIRECTED 

ACCELERATION 


Figure I. Einstein Equivalence Principle 


The experiments Umt led to this concept of space and time go back in history to 
the Renaissance, when Galileo Galilei is alleged to have tested the universality of free 
fall by dropping an iron cannon ball and a wooden ball from the Leaning Tower and 
observing that they fell in the same manner. Although this particular test is probably 
i^ocryphal, It is true that Galileo performed experiments to observe the rate of fall 
of objects mede of different materials by using inclined planes to dilute the effects of 
gravity. He also made some observations on pendulums, but it remained for Isaac 
Newton to employ this technique in a more careftil and quantitative manner. 

At this point, it was recognized that there might be a conceptual difference 
between the acceleration due to gravity on a particle and the usual mechanical accelera- 
tion. The attribution of gravitational mass and inertial mass was given to the material 
body suen that m^a = F and m^g = F, and the question of the proportionality of m^ to m^ 
regardless of the composition of the mass resulted in experimental tests with pendulums 
of various types. In 1686, Isaac Newton published his monumental Mathematical 
Principles of Natural Philosophy - or Principla - and the very beginning of the English 
translation from Newtonian Latin by Motte in 1729 refers to a set of pendulum experi- 
ments "very accurately made;" the cover page and first page of this work are reproduced 
as Figures 2 and 3. 


These experiments were later repeated by Bessel with a reported accuracy 
of about 1 part in 1000. Then, in the 1890, Eotvos (1922) conducted a series of tests 
using a torsion pendulum of very long period to look for the effect of the earth's gravity 

(m effect) and the effect of centrifugal acceleration (m. effect) on various kinds of 

® -9 ^ 

masses. He reached a level of about 10 in the constancy of the ratio of m, to m . 

i g 

These data were undoubtedly a large influence on Einstein's thinking during the forma- 
tion of his theories in the early 1900s. This type of pendulum e^qjeriment was also 
repeated with far greater precision using the sun's gravity at the earth's surface in the 
1964 to 1971 time frame by R.H. Dicke (see Roll et al . . 1964) and by V. B, Braginsky 


and V.I. Panov (1971), whose experiments achieved an accuracy between 10 

10 - 12 . 


■11 


and 


The formal theoretical extension of these tests for matter as a generalization 
of the Einstein equivalence principle is still in progress, and various forms of the 
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Figure 2. Newton's Principla. 
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MATHEMATICAL PRINCIPLES 


OF 

NATURAL PHILOSOPHY* 


Dtf initio ns 

DEFINITION I 

Tfu quantity of matter it the meature of the tame, anting from itt ^entity 

and haH conjointly.* 

T hus ah of a double density, in a double space, is quadruple in quan- 
tity; in a triple space, sextuple in quantity. The tame thing it to be 
understood of snow, and fine dust or powders, that are condensed 
by compression or lique&ction, and of all bodies that are by any causes 
whatever differently condensed. I have no regard in thu place to a medium, 
if any such there is, that freely pervades the interstices between the parts 
of bodies. It is this quantity that I mean hereafter everywhere under the 
name of body or mass. And the same is known by the weight of each body, 
Ibr it is proportioRal to the weight, as I have found by experiments on pen- 
dulums, very accurately made, which shall be shown hereafter. 

DEFINITION II* 

The quantity of motion it the meature of the tame, anting from the 
pelocity and quantity of matter ron jointly. 

The motion of the whole is the srm of the motioru of all the parts; and 
therefore in a body double in quantity, with equal velocity, the motion is 
double; with twice the velocity, it is quadruple. 

I* A» w Si « .Nolt IP 1 l'Apprada.Nal(ii.] (• Apptadii. Nmcu.] 

CiD 


Figure 3. First page of Motte's translation. 
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^ equivalence principle have been devised. The so-called "weak" equivalence principle 

applies only to conclusions obtained from tests for matter, the "strong" principle being 
^ Einstein's. By invoking other physical laws, such as the conservation of energy, it 

[ can be argued that connections exist between the two. However, like all the laws of 

[ physics, the principle of equivalence has its ultimate Justification in experimental test. 

( 

The generul term "Equivalence Principle" is often referred to in terms ot the type 
of data that support it. The Eotvus-type eiqieriments eiq)resslng the universality of 
free-fall trajectories are referred to as supporting the Weak Equivalence Princ4>le 
(WEP). The redshift phenomenon tests the Einstein Equivalence Principle (EEP). 
Schiffin I960 argued that the WEP implies the EEP for any self-consistent theory 
of gravity. Lightman and Lee (1973) were able to prove Schiff's conjecture for 
electromagnetically interacting systems in a static spherically symmetric gravitational 
field. Ni (1977) later showed that by loosening their assumptions, violations of Schiff's 
conjecture could be found for spinning test bodies. 

By removing the implicit assumption that material bodies are governed by the 
same metric as electromagnetic fields, Ni (1979) has recently investigated what is 
proved by certain aspects of different types of Equivalence Principle experiments 
typified by those of Eotvos, Hughes (Hughes and Robinson, 1960) and Drever (1961), and 
the redshift experiment being reported here. 

Ni finds that the Eotvos and Hughes-Drever ejqieriments demonstrate the equality 

1 2 

of the space-space components of the two presumably different metrics at the 10~ 
level. To describe the equality of the space-time components this level of accuracy 

_4 

must be reduced by a factor v/c ~ 10 , where v is the earth's orbital velocity, to 
-8 

yield a test at the 10 level. In the case of the time-time component comparison the 

-12 2 2 -4 

10 accuracy level must be reduced by a factor v /c leading to a test at the 10 

level. The present, direct test of the Einstein Equivalence Principle is aimed at this 

level of accuracy. 
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3. THE GRAVITATIONAL REDSHIFT AS A TEST OF THE EQUIVALENCE PRINCIPLE 


The simplest demonstration for the predicted connection between the effect of 
gravitation on time ^redshift) and the equivalence principle can best be given by the 
following type of argument, which is illustrated in Figure 4. 


Let us assume that a radio transmitter is located at the east wall of a laboratory, 
with a receiver on the west wall a distance i away. The signals will take a time 
At = i/c to cross the laboratory, c being the velocity of light. Now let the laboratory 
be accelerated eastward by a bulldozer charging into the west wall. During the time 
interval At, the laboratory will be accelerated to a velocity aAt = ai/c, and the fre- 
quency of the received signal at the west wall 'will be doppler shifted by a fractional 

2 

amount: Af/i = Av/c = ai/c . The question is whether, if we tip the laboratory on 

2 

edge east end up, the fractional shift would be Af/f = g£/c ; that is, can we replace 
a by g? This ejqjressio.: takes on the usual form for expressing redshift if we write 
gi as the Newtonian potential A4> and thus obtain 


f 2 • 

c 


( 1 ) 


We can answer the above question if we write the metric e:q)ression for an interval 
of space-time ds using the free parameters p and y, as done in the static and isotrcpic 
expression given by Eddington (1922) and Robertson (1962), and also including in the 
redshift a free parameter a. We have 

2 

here, M/R s Gm/c r^, where G is the universal constant of gravitation, m is the 
mass of the spherically symmetric gravitating body, and r is the distance from its 
center. 
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CONSTANT ACCELERATION = a 



Figure 4. Illustration of Equivalence Principle. 
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If we connect two fixed regions of space-time with light signals, keeping to terms 
2 

of order c , we have 



( 3 ) 


Since we are interested in the departure of a from unity, we can write a = 1 ± £, where 
±€ is a measure of the invalidity of the principle of equivalence. We should note here 
that the very concept of the existence of a metric description of space time depends on 
a being equal to 1 and that the use of the symbol in equation (2) (which is a metrical 
description) is therefore tantological. Nevertheless we choose to retain it and defend 
its presence as a means for testing the self consistency of the theory. 


4. THE BASIC CONCEPT OF THE GRAVITATIONAL REDSHIFT SPACE-BORNE 
EXPERIMENT 

The redshift experiment using clocks in a spacecraft and on earth compared by 
electromagnetic signals can be described in simple terms as follows: ’’Take a perfect 
(or "proper") clock and position it at several locations having a different gravitational 
potential and measure its rate at each position by comparing it to a clock that is at a 
constant gfravitational potential." By plotting the relative frequency difference versus 
gravitational potential, we can test the extent the slope of this relationship departs 
from imity in equation (1). The operation of comparing is done by means of light 
signals, either by looking at the hands of the clock in the literal sense or by transmitting 
time infonnation (or phase), to give the rate of time variation (or frequency) at a 
particular gravitational potential measured in terms of the standard or fixed clock. 

Here, the clocks are assumed to be totally unaffected by gravity fields or by accelera- 
tion. They are assumed to be the proper clocks of theoretical relativity. 

At present, the best embodiment of a proper clock for the purpose of such a test 

is the atomic hydrogen maser, whose stability, described in terms of the Allan two- 

sample variance (Allan, 1966), is shown in Figure 5. This form of variance, often 

referred to as the Kolmogorov variance, expresses the l<r expectation value of the 

frequency difference between adjacent measurements each of duration, t, in a time- 

14 

ordered series of data. Since the maser develops stability better than 1 part in 10 

for time intervals beyond 100 sec, we need at least 100 sec to get data for comparison 

-14 

with a Icr accuracy of frequency at the 1 X 10 level. 

The next requirement is to place the clock at different gravitational potentials. 

In this test using earth's gravity, we predict the redshift versus radial distance from 
the earth's center as shown in Figure 6. The idea of making comparisons at different 
gravity potentials can be implemented by using a nearly straight up and down ballistic 
trajectory and measuring the redshift frequency continuously as the vehicle travels 
in space. 
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Figure 6. Redshift vs. r?.dial distance from earth center. 


The original idea of using an eccentric orbit would allow the experiment to be 
repeated many times, thus improving the statistical validity of the data; however, under 
conditions such as this case where the rocket's weight-carrying capabity is small and 
the payload operating time is limited, it is better to put the rocket's energy into obtain- 
ing gravitational potential. In other words, the strategy is to reach as hig^ an altitude 
and to obtain as many data as possible across a large difference in gravity potential. 
Figure 6 shows roughly the altitude and redshift for a range of payload weight between 
200 and 220 lb (91 to 100 kg) when propelled nearly vertically by the four-stage Scout D 
rocket system. 

At this point, it is probably more realistic to refer to the clock as an oscillator 
and turn the discussion in terms of frequency rather than time or phase. While the 
distinction is mostly one of semantics, for clarity in pictorial discussions, we will 
use frequency, which is the parameter we are comparing as a function of gravitational 
potential. Our experiment, then, is to test the observed frequency of the probe clock 
in tlie manner shown in Figure 7 and establish the nature of the parameter a = 1 ± e . 

-10 2 

Since the test is expected to span a range of 4 X 10 in A<)>/c and our oscillators 

-14 

have stability greater than I X 10 throughout the ejq^ected time interval of the mission 

we can aim for a maximum sensitivity in the determination of the value of a of 

1 X 10”^V4 X lO”^^ or 25 parts per million. This excludes the question of systematic 

errors inherent to the ejqjeriment hardware, signal propagation and of our knowledge 

of the trajectory upon which our predictions will depend. Because the probe's trajectory 

is ballistic with .'in apogee altitude of about 10, 000 km, we will encounter large first- 

order doppler shifts in the micxx>wave signals connecting the space and earth station 

clocks. Ejqjressed as fractional frequency changes these will be on the order of 
-5 

1 to 2 X 10 , some nine orders of magnitude greater than the stability of the oscillators. 

“9 

The first-order doppler effect must be accounted for at the 10 level if we are to 
achieve the expected accuracy in the redshift measurement. 

In our considerations the possible error sources and the generation of error 
budgets to guide the design of the e;q)eriment, we insist that the errors in the oscillators 
should predominate and that the oscillators and their associated signal processing 
equipment will perform within or below the shaded band in Figure 5. It is obvious 


DATA REDUCTION 



Figure 7. Illustration of test for linearity of a. 
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that frequency contributions to error at the level of a few parts in 10^^ can mount up, 
even when considered to be statistically independent and combined in a root-sum- 
squares manner. For this reason, all frequency errors at the 2 X 10~^^ level or 
greater in our error budgets are considered to be significant. 
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5. THE DOPPLER CANCELLATION SYSTEM 


Next to the clocks, the most crucial part of the eiq)eriment is the system that 
removes the first-order doppler effects so that the relativistic data are revealed for 
measurement. This system is shown in very simple form in Figure 8, where a three- 
link microwave system connects the probe and earth oscillators. One link transmits 
the probe oscillator signals to the earth, where the comparison is made. This com- 
parison is shown as the beat signal from the microwave mixer Ml in Figure 8 and 
contains the first-order doppler and relativity information. A second microwave sys- 
tem attached to the probe receives a maser-contiolled uplink signal from earth and 
phase coherently retransmits it back to earth via a transponder. The returned signal 
is compared with the uplink signal and mixer M2 in Figure 8; the resulting beat 
frequency contains twice the first-order doppler-frequency shift, which, when divided 
in frequency by 2, is subtracted from the probe-maser signal at mixer M3, thus 
removing the first-order doppler frequency from the output frequency. 

The development of this system began with the original suggestion made in 1960 
by Badessa et al . Since we will require continuous operation of three microwave links 
each will require its own frequency to prevent self Jamming and regeneration. Figure 
9 conceptually describes the use of rational frequency multipliers to generate the 
separate frequencies required. 

It was pointed out to us by Prof. R. V. Pound (1973) that the use of three separate 
carrier frequencies in the S-Band region, spaced apart in frequency by the transponder 
input-output frequency ratio imposed by our use of the existing Unified S-Band system 
transponder, would produce substantial ionospheric doppler shifts that result from 
the rapid changes in columnar electron density in the signal paths. The system, shown 
conceptually in Figure 8 will cancel nondispersive propagation effects such as the 
refraction from the earth's troposphere; however to assess its behavior in a dispersive 
medium such as the ionosphere we must look at it in more detail. 

To bring the prqjerties of the doppler and ionosphere cancelling system together, 
we can discuss the total doppler effect Af^f in terms of the rate of change of the 
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Figure 8 


Doppler cancelling system concept. 



















"optical path" of the electromagnetic energy through the medium having a refractive 
index n. The atmospheric and ionospheric dc^pler shift resulting from propagation 
along a path P is as follows: 

p 

Here, n(t) has three components corresponding to a vacuum (n » 1), atmospheric 
refraction n , and ionospheric refraction n.. Of these, only the ionospheric term 
depends on frequency for signals near 2 GHz. Consequently except for ionospheric 
dispersion, all effects between the spacecraft and ground stations are cancelled pro- 
vided that the antenna system treats the signals at all three frequencies in the same 
manner. 

This point is crucial to the success of the technique, for it is the antenna that 
represents the trajectory of the payload, and its position is considered to be the point 
in space that interacts with the iqjcoming electromagnetic signal from the ground- 
station transmitter and simultaneously phase coherently retransmits this signal while 
simultaneously transmitting the clock downlink signal to the earth. 

For the moment, wp can consider that our system will strictly account for the 
probe motion as represented by its antenna, and the cancellation is complete if n = I . 
The nondispersive refraction from the atmosphere is cancelled exactly only if the 
path integral P of the upward signal from the ground station is exactly the same as 
that of the downlink signals from the clock and the transponder. These effects and 
those of the earth station motion during the light time of the signals will be considered 
later. 

To make use of existing Unified S-Band (USB) transponder designs, we adopted a 
turn-around frequency ratio of 240/2!' i lor the transponder uplink and downlink. The 
resulting s^aration of these frequencies requires a technique to cope with ionospheric 
dispersion. 
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Phase shifts that result from propagation through the ionosphere can be evaluated 
by tracing Uie signals shown in Figure 9. The refractive component of the ionosphere 
nj is a function of position { and time 1, and depends on electron density and magnetic 
field as follows: 


njU,t) = 



( 6 ) 


1 P 6* 

Here, f is the propagation frequency and the plasma frequency ~ ^ ^ depends 

on the number density of electrons p({,t), the electronic charge to mass ratio e/m, 
and the permittivity of free space €«. The frequency f is the electron gyromagnetic 

V m ^ 

frequency 2.8 mHz per oersted of magnetic field. The ratio fjj^/f is less than 1 X lO”** 
even for maximum terrestrial magnetic fields and this term can be ignored. Under 
these conditions following Tucker and Fannin (1968) we identify the do: it term of 

a power series e^qjansion in frequency f as follows: 


and 




ni({,t) = (l - 


4tt2 ‘o 


8 it^ ‘o f^ 




t) ds 
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(6) 


( 7 ) 


In the system shown in Figure 9 we can trace the frequency shifts resulting from 

■mJ^ 

changes in the columnar electron density p(5,t) and obtain the ionospheric frequency 
error at the input to mixer M3 
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Here, is the hydrogen maser output frequency (Hellwig et al . . 1970) 1,420,405, 751.68 
Hz and M/N, P/Q, R/S are numerical ratios used to synthesize the frequencies in the 
three microwave links. Similarly, for the spacecraft oscillator downlink we find an 
ionospheric frequency shift given at the other iiq>ut to mixer MS as follows: 


Af 


1 1 way 
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The mixer MS produces a net shift given by the difference of these two shifts: 


Af = Af. 
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(10) 


The ratios P/Q, M/N, and R/S can be chosen to make the last term equal to zero 
and, to the extent that the approximations made in equation (7) hold true, including 
the effects of ray path separation owing to dispersion and elapsed time between upgoing 
and downgoing rays, the effect of the ionosphere can be cancelled. 

Since our system made use of already existing Unified S-Band equipment, the 
trarisponder ratio, M/N, as mentioned earlier was 240/221. This predetermined ratio 
and the usable frequency range of the USB system led us to choose P/Q = 76/49 and 
R/S = 82/55. 

The residual error given by 

otj 2 \ n/r2 / ^2 

1 way 2R \ M / Q 
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when the ratios are chosen as above is 




Af„ 


aT 


= 2.5 X 10 


-5 


1 way 


For typical ionospheric conditions and with the trajectory flown we find that 
Af, Ji is, at worst, about so that Af /f is, at worst, 2.5X10 

The ionospheric frequency signature (Figure 10) shows the uncancelled behavior 
predicted from measurements of ionospheric columnar electron density made during 
the experiment. 


The extent to which the approximations made in equation (7) and the effect of ray 

bending and retardation (nonsimultaneous ray palhs) through the troposphere and 

ionosphere have been studied by Mr. C. Baugher at the Marshall Space Flight Center 

using more exact ray-tracing methods ( 1980). These values of the ionospheric 

profile were measured at 1145 and 1330 GMT, June 18, 1978. These calculations 

-■> -3 

were made to estimate the f and f terms not included in the expansion given 

-2 

in equation (7). Here the f term is identified with Faraday rotation of the plane 

-3 

of polarization; this is a nonreciproeal effect in the propagation direction. The f 
term is chiefly associated with the differential bending of the ray paths that result 
from ionospheric dispersion. The system, in fact, cancels the electromagnetic effects 
of ray paths that are not exactly coincident in time but are separated by a time inter- 
val of r/c sec. 


Taking these effects into consideration, the ray-tracing computations, including 

time delay, between uplink transmission and downlink reception, show a significant 

14 

effect Af/f > 0, 5 X 10 only during the last 2 min before impact. As will be seen later, 
the data taken after the uplink transmitter interruption, which occurred 12 min before 
impact, were not considered reliable and have not been used in the present reports 
of the redshift results. 


27 


^ LINEAR SCALE 


10 

9 


-ex 10' 



10 


r8 


10 


10 


>-io 


10 


.-II 


10 


-12 


10 


13 


10 


14 


10 


-15 


Figure 10. Uncancelled ionospheric frequency shifts for trajectory flown June 18, 1976. 
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Variations in the troposphere and ionosphere that occur rapidly compared to the 

light time (which is always less than the value at apogee of 0. 03 sec) will produce 

phase variations that are not cancelled by this system. These will occur only at 

Fourier frequencies higher than 33 Hz and since our data averaging intervals to obtain 

-14 

frequency stability at the 10 level are greater than 100 sec, these short-term phase 
noise fluctuations will not affect the experiment. 
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o. PHASE STABIUTY CONSIDERATIONS IN THE PROBE AND GROUND STATION 
ANTENNA SYSTEMS 

A very critical factor in the frequency stability of the system shown in Figure 9 is 

the phase stability of the electronic systems, cables, and propagation paths. This is 

-14 

vividly brought into focus when we consider that a fractional frequency shift of 1 X 10 

-2 

at 2.2 GHz results from a change of only 1.3 X 10 rad in 100 sec or a change in cable 
length of 0.03 cm in the same time interval. Questions relating to the dispersive and 
nondispersive aspects of the propagation have been discussed and we have seen that 
all perturbations resulting from the relative motion of the antennas are cancelled within 
the light-time interval r/c. However, this cancellation does not apply to the cables 
connecting the ground station antenna with the transmitters and receiver parts of 
Unified S-Band system (USB) located some distance away. In the normal USB system 
the two receivers share a common low-noise "front end" and the signals are carried 
via common cables to the phase-locked tracking filter receivers. The uplink trans- 
mitter signals are sent to the antenna along a separate cable. Because tte doppler 
cancelling system subtracts one-half the two-way doppler signal from the one-way, 
only one-half of the phase fluctuations in the receiver cable are cancelled and none of 
those variations of the transmitter cable are accounted for. In addition to this totally 
uncancelled aspect of the transmitter cable there existed the very large uncontrolled 
phase variations in the signal path from the transmitter exciter, through the many 
stages of frequency multiplication and on to the high power amplifiers. This problem 
was solved by using a separate, specially controlled, cable to sample the transmitted 
power near the feed point of the antenna and bring this signal back to the control 
station where the hydrogen maser master oscillators were located. The phase of this 
signal was compared with the phase of a maser-generated signal, and the offset 
signal, proportional to the phase difference, was used to control the phase of the 
transmitter master oscillator. In this way the fluctuation of cable length, power 
amplifier and frequency multiplier stages were all placed within the phase-locked loop 
and the output signal. 
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In the case of the spacecraft antenna, transponder, and downlink transmitter, 
the situation was very different. Here cable dirtances were short and the transponder 
and transmitter operated their three separate links through a triplexer that routed 
incoming signals to the receiver and sent both outgoing signals to the antenna. 

The chief concern in the probe was the phase stability of the transponder and 
downlink transmitter during the large temperature changes expected when the space- 
craft was launched into space. Since no time could be allowed to restabilize these 
components while in space their temperature control systems were required to cope 
with unusually large perturbations. The usual USB transponder design incorporated 
two separate units within a single enclosure. For the present application, the USB 
transponder was modified by replacing one ti^nsponder section with the frequency 
translator system, shown as the P/Q block in Figure 9. The complete transponder and 
translator package were temperature-controlled by a liquid ammonia boiler, which kept 
the temperature within a range of 0.25°C during the part of the mission when useful data 
were acquired. Our strategy for the one-shot mission was to stabilize the device 
thermally for several days before the flight and then to launch it into space with the same 
thermal-control system used on the ground. The thermal control system, specially 
developed at the Marshall Space Flight Center for this e:q)eriment, consisted of an 
ammonia boiler that contains enough liquid to last the duration of the mission. The 
boil-off temperature was controlled by gas pressure, which was maintained by a 
venting orifice. Remarkably close control of temperature was achieved in both the 
transponder and the translator; these units were located in the same insulated con- 
tainer housing the USB transponder, which was in thermal contact with the boiler. 
Performance data and the calibrations for phase with temperature are discussed in 
Section 14. 

Phase variations resulting from the antenna radiation patterns are also important. 

In the case of the ground station antenna, which is steered toward the probe, the 
signal propagation directions is very constant with respect to the axis of the antenna. 

In the case of the probe, which did not have a steerable antenna, the signal phase 
depended on the relationship between the signal propagation vector and the direction 
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of the antenna dipole; the phase shifts also depended on the frequency of the signal in 
each of the three links. 

Since the probe was spin-stabilized, the axis of rotation was chosen as the center 
of symmetry of the dipole antenna system. The antenna, Including its ground plane, 
in fact encompasses the entire exterior conducting surfaces of the vehicle. The space 
vehicle was therefore completely enclosed in an electrically conductive shroud that 
was made S3rmmetric in rotation about the spin axis. 

The dipole antenna that operated the three links was located on axis in a 10* 
conical ground plane. Under these conditions, rotation of the payload will cause no 
change in signal phase. Because the location of the phase center of the antenna is 
slightly different for each frequency some change of phase with aspect angle is 
expected. Careful atenna design was required to provide a low level of phase varia- 
tion. 


Figure 11 shows the predicted combined effect from the three signal links when 
processed by the system shown in Figure 9. 
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7. PREDICTIONS MADE FROM RELATIVITY THEORY 


This section discusses the behavior of the output frequency from the system shown 
in Figure 9 according to the predictions made from the general theory of relativity. 

The derivation follows the notation and method used by Kleppner, Vessot, and Ramsey 
(1970) to obtain the relativistic prediction of the one-way and two-way doppler shifts, 
which are then combined by subtracting one-half the two-way signal from the one-way 
signal to obtain a doppler-cancelled signal containing only relativistic effects. 


We choose an inertial earth-centered coordinate frame as shown in Figure 12 and 
apply the doppler equation of general relativity to an oscillator having a proper fre- 
quency fg at a position moving with velocity measured by receiving a signal of 
frequency f ' at a position Tg moving with velocity Pg. The unit vector describing the 
inertial propagation direction from the position to the rg position is Cgg. Keeping 
terms of order c"^, we obtain f', the frequency at the ground station received from 
the probe: 



1 2 “^2 

1/2 

" ^ * 

II 

0 

1 + (2(j>2/cl - p“ 


^ * ‘23 

1 +(2<t>g/c^) 


■ ^2 * ‘ 23 . 


( 11 ) 


where (|)2 and (|)g are the gravitational potentials of the probe and the ground station, 
respectively. If we expand this expression binomially, the term in the first bracket 
is recognized as containing 1 + (( 4>2 - 4 * 3 )/ c^] - KI^I - |^|)/ 2 ], which is the com- 
bination of the redshift and second-order doppler effects. The second bracket in 
equation ( 11 ) contains the first-order doppler effects, which are some 10 times larger. 

In a similar manner, we can successively apply the doppler equations of general 
relativity to the transponder signals: 
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( 12 ) 



l4(2<t>/0-p‘ 

1 + ( 2 <|> 3 / c 2 ) -tl 




V ‘^1 ’ ‘l2/ 


Here we observe that the subscripts I and 3 refer to the ground station at the time of 
transmitting and receiving the signal. During this time (which equals 2r/c, where r 
is the path distance to the probe), the earth has rotated and ^3 has changed direction 
but not magnitude. The earth's gravity has not changed, so = 4 > 3 « 1'he term in the 
first bracket thus cancels to 1 , leaving only the first-order doppler expression in the 
two remaining brackets, which, if we £q>proximate by putting = P^and Cgg = “* 12 ’ 
as expected, converges to twice the doppler frequency: [1 + 2 (p 2 - Pg) • C 23 J, given 
by f" - f^. [Note: This last approximation is given only to illustrate magnitudes and 
will not be applied in the following derivations.] 


The system in Figure 9 combines the two downlink signals to obtain at mixer M3 




and 


f f" + f 

Af _ *0 *0 ^0 
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(13) 


Substituting equations (11) and ( 12 ) into equation (13) and setting = <j>g and |p^ | = |Pg|, 
we obtain, after considerable manipulation, 
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where, in the second-order terms, we define 


6 


e 




and 
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This e;q}ression can be put in the following form, which we use to describe the relative 
motion between the earth and the probe: 


At *^e 1*^2 " At— — 

T ~ 2 ^ ‘23 * ‘"S • 

c 


(15) 


We note, however, that and are not the same, since the earth, whose 
angular rotation rate is SI, has rotated an angle QAt, where At is the go-retum time 
of the transponded signal. The direction of p^ has changed, and the corresponding 
acceleration representing this velocity change is the familiar centripetal acceleration 
a in the radial direction perpendicular to the earth's axis of rotation. 

The last term in equation (15) expresses the fact that the earth moved during the 
light time At of the two-way signal and that the doppler cancellation system, which 
cannot react faster than the speed of light, must account for the change in velocity of 
the earth's rotation during the light time of the signal. This is an exact replica of the 
situation shown in Figure 4, which was used to Illustrate the redshift. Here the 
acceleration is the result of the rotating earth and, assuming the coordinate frame 
we are using, would vanish in our experiment only if we chose to locate the tracking 
station at one of the earth's poles. The residual effect is the first-order doppler shift 
of a receiving station, whose velocity changes during the time At/2 (one-way owing to 
an acceleration a^). The component of this velocity change in the propagation direc- 
tion is represented by the dot product *23 • Qg, where (^2 vector describ- 

ing the direction of propagation from the probe to the earth station. Figure 13 shows 
the magnitude of this term as a function of time for a trajectory similar to the one 
actually flown. 

The major relativistic effects are contained in the first two terms of equation (15). 

The first is the redshift effect we seek to test, while the second is the second-order 

doppler effect of special relativity due to the relative velocity of the earth and the probe 

oscillators. It is interesting to note that, because this e?q)eriment is essentially one 

2 2 

involving a freely falling object, the contribution of the v /c second-order doppler 
effect (where v is the velocity) is equal to the redshift. This can readily be seen if we 
consider the conversion of potential energy <|) to kinetic energy over the trajectory 
between points A and B: 


Af/f IN UNITS OF 10 









( 16 ) 


B 



This is a manifestation of the virial theorem of classical mechanics, where the total 
energy is equal to the average value of the kinetic energy. 

The result of combining the redshift and the second-order doppler effects is shown 
in Figure 14, which are very nearly equal as a result of the free fall of the probe. We 
observe there the effect of the negative sign of the second-order doppler effect at the 
ends of the trajectory, where the relative velocity between the earth and the probe is 
greatest, and the effect in the positive direction near apogee, where the redshift term 
dominates. 


The objective of this experiment is to see how A^f varies with gravritational 
potential. We must therefore account for the second-order dq>pler and earth motion 
effects shown as the second and third term of equation (15) in addition to predicting 
the gravitational potential variation shown in the first term. Position, velocity, and 
direction information are required as a function of time at an accuracy commensurate 
with the clock stability. In addition to this we must know the gravitatio.ial potential 
including possible effects of lunar and solar gravity as a function of the probe's posi- 
tion. 


Knowledge of the time dependence of the gec^otential at the probe, along with 
position and velocity data, is used to generate a prediction of the right-hand side of 
equation (15). The comparison of the prediction with the actual data Is, in fact, 
the test we are making. 

The fact that the microwave signals traverse large distances over which substantial 
changes in gravity field occur means that this comparison of oscillators is a measure- 
ment of very nonlocal effects. Furthermore it relies on the assumption of the isotropy 
of the one-way velocity of light. To predict the behavior of the data from this rather 
complex system we will assume that the oscillators, electromagnetic signals, velocities, 
and gravitational fields are all interacting in the self-consistent marmer described by 
Goieral Relativity. In short, this esqjeriment tests the internal self-consistency of 
these aspects of General Relativity. 
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8. PROBE TRACKING REQUIREMENTS 


The experimental test of the gravitational redshift depends on our knowledge of 
all the terms in equation (15) at all points in the trajectory. The data-reduction pro- 
cess we use will compare the measured value of Af/f with the predicted value of At^f 
obtained from our knowledge of the geopotential, the relative velocity, and the value 
of the component of the acceleration of the ground station due to the earth's rotation 
in the direction of the line of sight to the probe. We will presume that the second- 
order doppler and earth-acceleration terms of the special theory of relativity are , 
in principle, correct and will assign all the discrepancies in the data to the gravita- 
tional effect ((|)p - <|>g)/c^. 

In keeping with the constraints on accuracy set forth earlier we require that the 

overall error in the trajectory determination must be such that the error in the red- 

shift prediction, including the second-order doppler and earth rotation terms of equation 

-14 

(15), is less than 1X10 in Af/f. To achieve this, we set the error in the position 
determination, which, in turn, governs the gravitational potential prediction as follows: 



-15 

For an error contribution of ±2 X 10 near apogee (altitude 10,000 km) we must 
determine position in the radial direction within ± 100 m. For the same error near 
earth, the radial distance must be ± 16 m. 


The error contribution from uncertainties in the second-order doppler are as 
follows: 
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For the same level of error taken as ±2 X 10 we must determine the velocity 
within 6 cm/sec near apogee and 1 cm/sec near earth. The uncertainty in the angular 
position of the centrifugal acceleration in relation to the propagation vector is consis- 
tent with the position requirements given above. 

Under the above position and velocity restrictions on the tracking accuracy, the 

root-sum-square contribution of tracking-related parameters to the experiment's 

15 

accuracy is about 3. 5 parts in 10 . 

The question of position and velocity tracking accuracy is crucial to the accuracy 
of the mission. At the beginning of the program, two separate studies were made to 
ascertain the feasibility of the experiment in this regard (Townsend and Alford, 1972; 
Ryan, 1973). The recommendations were to use range-rate data exclusively and to 
take these data from at least three tracking stations equipped either with atomic 
hydrogen masers or improved cesium-beam frequency standards. 

The USB system had four stations that can provide useful data; their station loca- 
tions are as follows: 



Station 

Latitude 

Longitude 

1. 

Merritt Island, Florida 

28'30'30" 

279“18'24" 

2. 

Bermuda 

32 21 5 

295 20 32 

3. 

NTTF 

38 59 55 

283 9 25 

4. 

Ascension Island 

-7 57 17 

345 40 23 


Of these, stations 1, 2, and 3 were equipped with hydrogen masers; the masers at 
Bermuda and the National Test and Training Facility (NTTF) were built by the Goddard 
Space Flight Center (GSFC). The Merritt Island Laimch and Acquisition (MILA) station 
was the primary data-acquisition site, recovering the relativity data with the system 
shown in Figure 9. Two SAO-constructed VLG-lOA masers were located at MILA to 
serve several purposes: to provide redundancy in the event of failure, to furnish an 
internal self-checking capability throughout the actual mission, and to enable preflight 
simulation tests to be done before the mission. 
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Doppler tracking was done at Merritt Island using both the two-way and the one-way 
methods. The two-way doppler data were recovered by measuring the doppler beat 
frequency between the uplink and the translated down-link signals. The one-way 
doppler data, which also included the relativistic frequency shifts, were obtained by 
comparing the clock signal received from the probe with the station oscillator. At 
MILA both t 3 q>es of data were logged in ten.is of actual counts of doppler-frequency 
zero crossings taken over 0.2-sec intervals. The one-way data were taken at the 
even tenths of a second, the two-way at the odd tenths. 

Doppler tracking at the other stations was done by logging three-way doppler 
cycles, the latter resulting from the behavior of the frequency of the signal that 
originates in the uplink transmitter at Merritt Island and is transponded to the several 
ground stations. 

The operations plan allowed for a minimum of three stations in addition to Merritt 
Island to be ready at the time of launch; the Merritt Island, Bermuda, and NTTF 
stations were considered as the primary tracking stations, with Ascension Island 
serving as back-up station Merritt Island had the additional responsibility of being 
the redshift data-acquisition station. All stations recovered telemetry data from the 
probe. 
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9. REDSHIFT DATA ACQUISITION 


The primary redshift data from the last mixer shown schematically in Figure 9 
were recorded in digital form by measuring the voltage of the sinusoidal signal at 14-bit 
accuracy (0. 006%) every 0. 01 sec. Output data in the form of both sine and cosine 
voltages were acquired so as to resolve ambiguity of phase reversal, which, as men- 
tioned before, occurred at two points in the mission. 

The instability that these discontinuous and granular measurements produced in 
the data are shown in Figure 15 and are well below the instabilities of the maser 
oscillators being compared. A redundant data-acquisition system operating completely 
independently of the primary system was also included at the MILA site to act as a 
backup in the event the primary data system malfunctioned. This system sampled the 
data at the same rate but at 12-bit granularity. 

The MILA station data were acquired and placed into the 1-frame/ sec format 

depicted in Figure 16. To relate the redshift beat-frequency data to the predicted 

-15 

redshift data with a maximum error of 2 X 10 , the time S)Tichronization between 

the tracking stations and the MILA redshift data-acquisition station must be maintained 

-4 -5 

within 10 sec. Since a synchronization accuracy at the 10 -sec level is routinely 

available to satisfy the tracking requirements; the e^qjerimental error contribution 

from this source is negligible. 

Timing for the 0. 01 -sec interval between 14-bit samples of the redshift data is 
done by freezing the voltage in a sample-and-hold amplifier on command from the 
100-pulse/ sec signal obtained from the station clock, which is monitored from Loran C 
signals controlled by the United States Naval Observatory. These signals represent 
the UTl time scale and provide the data required to determine the angular position of 
the rotating earth with respect to the usual astronomical coordinate axes. This 
information is necessary, since equation (15) is given in terms of the relative position 
of the earth and probe, and the probe trajectory, in terms of a coordinate system based 
on an earth-centered inertial frame. In representing the predicted values for equation 
(15), care must be taken to place the earth and the probe in the same coordinate. 
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Figure 15. Noise from recording system expressed as Allan variances 
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10. THE DATA-REDUCTION METHOD FOR THE REDSHIFT EXPERIMENT 


The actual behavior of the output frequency obtained from the doppler cancelling 
system is compared with the theory according to the flow plan shown in Figure 17. 

The left-hand arrows symbolize the three input data sources, which come from the 
format described earlier. These sources are 1) the sine and cosine output from the 
doppler cancellation system, 2) the probe trajectory as recovered from the doppler 
tracking data, and 3) the telemetry data describing the environmental conditions at the 
probe and the internal operating conditions in the probe maser, the transponder, and 
the translator. 

The overall program predicts the behavior of the redshlft beat signal from the 
trajectory by using the general theory of relativity as the working hypothesis. Further 
predictions relating to the previously calibrated phase and frequency sensitivities are 
added to the redshift prediction, including the fixed offset measured by direct com- 
parison of the probe and ground masers. These data are e^qjressed in terms of phase 
by integrating frequency data from an arbitrary epoch and encoding them magnetically 
to produce a data prediction tape. Because the trajectory information and phase 
correction data are relatively noise-free these data show little or no short-term 
variations. 

The sine and cosine data, are left unfiltered so as to recover the noise and what- 
ever rapid systematic effects resulted from the probe’s rotation and nutration. Values 
of redshift beat data are recorded in terms of phase versus time, starting at the same 
arbitrary epoch that was chosen for the prediction tape. As before, these data are 
recorded on magnetic tape to provide a redshift data tape. 

A 300-sec stretch of data centei'ed at apogee is used to determine the fixed offset 
frequency of the probe maser that resulted from the trauma of launch. This fixed off- 
set frequency is removed from the whole data set so that only the variation of the fre- 
quency from apogee is considered. In the results we find that the offset was very 
small, -3. 1 X 10 
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This offset information was intended for use in estimating the cavity frequency 
shift, which, if appreciable flux variation had occurred, would have caused systematic 
output frequency shifts. This correction was not required. 

In the data- reduction process all comparisons of prediction with actual data are 
made in terms of phase, thus avoiding the problems that could result if we attempted 
to determine, as a function of time, the frequency information in the data by filtering 
and differentiating. Filtering and differentiating that is required is done only to repre- 
sent the departure in the redshift data from the prediction for the statistical regression 
process, and since the frequency variations are small, no problems are anticipated. 

The geopotential at apogee is taken as the zero potential for our comparison of 

the redshift, A4>/c , with Af/f. In the data-flow diagram, the linear regression of the 

residual frequency, obtained by subtracting the prediction data from the redshift data 
2 

versus aVc , should have a slope of zero if the equivalence principle holds. Since 
there are random noise processes well as possible residuals, the data are repre- 
sented statistically as a slcpe c with a lo- uncertainty in the slope, ±6. Other fre- 
quency uncertainties resulting from the calibration and measurement errors determined 
during testing are computed along with the corrections applied to the prediction data; 
since they are uncorrelated, they are combined in a root-sum-square manner to repre- 
sent the l(T-accuracy capability of the experiment system and is designated as ±i. 

The combination of ±6 and ± i, again as root sum squares, comprises the combined 
Icr-accuracy capability of the experiment. 

In the development of the ejq)eriment plan, the accuracy capability was discussed 
in terms of a predicted error budget, which was used as a guideline for allocating 
errors among the various components of the experiment that result from random effects 
and systematic uncertainties. In the final analysis of the data, the random or noise 
effects can be compared to ±6 and the systematics with ±l to obtain an evaluation of 
the system performance. 

Additional means are available to evaluate the data and extract further systematic 
and statistical information. The residual phase data are treated as a time series and 


are processed to yield the Allan variance of the frequency comparison. This technique 

gives information about the spectral densities of the various noise components in the 

data, enabling us, to identify physical processes that can cause these disturbances. 

As shown in Figure 5, we see that the hydrogen-maser performance instabilities are 

caused chiefly by two processes: the behavior of the slc^e from 0.01- to 100-sec 

- 1/2 

averaging intervals and of the x ' slope for 100-sec averaging intervals and longer 

then an upturn resulting from systematic effects. These two slopes are characteristics 

of the additive white phase noise (x”\ the flicker-of-phase noise (x"^), and the white 
- 1/2 

frequency noise (x ) contribution originating in the hydrogen maser. 
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11. THE SPACE PROBE MASER OSCILLATOR 


Since the probe maser was specifically developed for the redshift probe experi- 
ment, and since its performance and our knowledge of its systematics under environ- 
mental stress are crucial to the outcome of the experiment, the construction of the 
maser is described here in some detail, together with a discussion of the major per- 
turbing environmental effects. Section 14 contains further details on the actual cali- 
bration data used to correct the redshift data for systematic variations owing to the 
^ace environment. 

From the data given for the ground-maser performance shown in Figure 5, we 

see that the contributions to errors from this source are on the order of 8 X lO*^^ and 
-15 

3 X 10 for 100-sec and 1000-sec averaging times and that, following the original 
precepts discussed in Section 4, the maser errors are taken as the chief limitation to 
accuracy. Since there has been no previous eTqperience, we e;q>ect that the error 
contributions in the space maser will dominate. The Identification of these errors and 
the removal of their systematics through tests and calibrations have constituted a 
significant portion of our effort in the program. These tests and calibrations were 
performed as part of the qualification test program for the actual flight hardware. 

This is the ’’protoflight" concept under which the program was governed. 

It is clear that there is a dilemna inherent to this conc^t. How far should a new 
device be tested to gain assurance of its mechanical, vacuum, and structural integrity 
without its flight- worthiness being lost? Many questions of this type were answered 
by testing an earlier engineering model of the flight maser, which v<as designated as 
the engineering verification unit (EVXJ). 

The list of tests of the EW and the later tests of the flight article will not be 
spelled out here. Detailed procedures for all these tests are on record. Many of the 
tests were substantially complicated by the requirement for accurate simulation of 
launch conditions for calibration purposes. Also, simulations of vehicle spinup and 
magnetic profiles during thermal -vacuum tests because of magnetic hysteresis effects 
had to be made in the proper sequence. 
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However, it must be emphasized that the success of the experiment resulted almost 
entirely from this massive test program, which served simultaneously both to qualify 
the probe maser and the other probe components and to calibrate the systematic varia- 
tions of phase (or frequency) resulting from enviroi omental variations from launch to 
the end of the mission. 

In addition to all the scheduled aspects of the test program, many unexpected 
problems arose, which had to be solved by making engineering changes, often on the 
spot, without perturbing the planned test schedule. It is clear that engineering 
develt^ment tests were frequently being made coincidental’ ' with the qualification and 
calibration tests. Required for all this was substantial flexibility in the test philosophy 
as well as good engineering practice. During this testing phase the problems that arose 
were solved in such a manner that the validity of a qualification result or a parameter 
for a calibration was not compromised by a subsequent engineering modification. Under 
the leadership of the MSFC program manager, the SAO and MSFC science, engineering, 
and test groups worked very effectively together to overcome many problems that often 
seemed impossible to solve within the schedule leading to the launch. 

11.1 Probe-Maser Design Requirements 

From the dynamics of a body falling nearly vertically in the gravity field of the 
earth and from the known behavior of atomic hydrogen masers, we can establish a 
very rough optimum situation, based on the following considerations: 

1. Allow enough time aloft to stabilize whatever launch-induced perturbations 
occur and minimize the effect of such thermal and mechanical perturbations so that 
the maximum possible stable operating time is available. 

2. Maintain communication between the payload and the groimd station throughout 
the mission; concentrate on obtaining data from apogee down to altitudes as low as 
possible at the end of the mission. 

Owing to the R“^ gravity potential acting on a freely falling body, the time avail- 
able for measurements near impact is a relatively insensitive function of the time aloft. 


56 


and the strategy therefore is to obtain a minimum period of about 2 hours, in which 
apogee redshift values greater than 4 parts in 10^^ are attained. The upward branch 
of the trajectory \vould also contain some data; however, because of the propulsion 
phase and payload stabilizing time, these data will be of less value than the data from 
the later, downward branch. 

Obviously, no time is available for the more conventional spacecraft outgassing 
and thermal stabilizing. Furthermore, since this is a one-shot mission, both the 
space and the ground equipment must (^erate without interruption and with the required 
stability throughout the mission. A very stable operating condition before launch and 
thermal and mechanical stability through the experiment are necessary. This plus the 
weight limitations resulting from the vehicle constraints placed rather stringent design 
requirements on the experiment. 

The first rrf;. irement is that of survival through the vibration, shock, acoustical 
pressures, and decompression of launch by the four-stage, solid-fueled Scout system 
into a zero g and spinning condition while in space. Of the total payload equipment, 
the hydrogen maser and the newly designed ammonia cooling system for the transponder 
are the onl}' types of items that had not previously been flown in space. 

At the outset, we allowed about 90 lb for the maser and adhered as closely as 
possible to this limit. The requirements for its survival and operation in space were 
as follows: 

1 . The maser-frequency fixed-frequency offset resulting from the trauma of 

12 

launch should be small, less than 5 parts in 10 , stabilizing within about 10 min to an 

14 

overall stability of 1 part in 10 for the remainder of the mission. 

2. Thermal, zero g, spin, and pressure effects from the transition into space 

-14 

causing possible longer term instability must be less than 1 X lO either by calibrat- 
ing for various environmental factors during preflight simulation tests or by engineer- 
ing the maser to cope with its immediate environment at this level of stability. 

3. The maser should be built so as to operate continuously for about 9 months 
to allow time for its qualification testing and calibration and to permit it to c^erate, 
without alteration, as a piece of flight hardware; this includes several weeks of pre- 
flight stabilization. 
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11,2 Maser Fiequency Perturbations 


The principal systematic frequency variation of hydrogen masers is described 
by the following ejqjression: 


f 


J_ 

f 



(f - f J + 2750 
' c 0' 


(17) 


The first term is the cavity-resonance mistuning, or "pulling, " effect, and the second 
is the second-order magnetic-field dependence of the atomic hydrogen hyporfine 
transition F = 1, Mj, = 0 — F = 0, M^, = 0. is the cavity resonator Q, and is 
the Q of the atomic transition, which depends on the geometry of the hydrogen-maser 
storage bulb, the quality of the wall coating, and the collision rate of the atoms among 
themselves (spin-exchange processes). This last process is a fimction of beam input 
flux, which, in turn, can be represented in terms of the maser output power level. 


Two aspects of the cavity-resonance shift especially concern us: 1) the variation 

in Af during the mission, and 2) the average magnitude of Af as a result of the com- 
c c 

bination of shake, shock, and zero gravity that occurs from earth-bound conditions 
through to the free-fall condition after powered flight ceases. The effect of cavity- 
resonance variations is obvious in equation (17). However, if there ir also a large 
fixed offset (f - f«), we are further subject to output frequency variations due to 
variations in resulting from changes in atomic hydrogen flux during the mission. 
These changes in flux result in changes in power level W, and the measured power 
level is used as a measure of the beam flux for calibration. 


In the probe maser, the cavity- resonance frequency is subject to many changes - 
from stress changes resulting from variations in the gas pressure P of the enclosure 
surrounding the cavity vacuum system, from thermal ejqjansion or contraction result- 
ing from changes in temperature T, owing to centrifugal forces resulting from varia- 
tions in the payload spin, and from changes in the cavity shape. 

The expression describing all the known perturbations to the output frequency of 
the maser is 


. ( 18 ) 


f 


fQ (M 
c I c 

, W 


dP + 


Of 
c 

9T 


dT + 


Of 
c 

an 


di2 


Q Af 
c c 


9U/QI 

aw 


dW + 


8B. 


dB, 


0 


In the magnetic calibration of tlie maser we observed that there was very little correla- 
tion between the second-order magnetic field sensitivity of the atomic hydrogen hjper- 

2 

fine separation pi'edicted from the expression 2750 H and the measured output 

frequency shift. At first we attributed this anomal^ to the cavity microwave isolator, 
which, bt'ing a forrite device operating in a magnetic field carefully adjusted to provide 
the proper resonance despite being shielded, was though to be affected by external 
fields. Subsequent developments and research on hydrogen masers disclosed that 
the effect was the result of changes in the state distribution of atoms among the 
h>perfine magnetic sublevels while in transit from the hexapole magnet to the bulb. 


Changes in population among the atomic hyperfine sublevels produce changes in 
the ou4>ut frequency of the maser as a result of interatomic spin-exchange collision 
processes duruig the storage lifetime of the atoms in the bulbs. Our solution to this 
problem was not only to calibrate the magnetic sensitivity but to simulate the whole 
magnetic history from prelaunch erection of the Scout rocket on through launch, 
spin up, and during the mission. 

The magnetic field behavior was measured by simulating the magnetic history 
of the payload before launch, during launch, and throughout tlie flight. A standard 
initial-condition state magnetization of the magnetic shields was ensured by a 
pregaussing operation performed before the vehicle-erection phase of both the actual 
flight and the magnetic simulation of the flight. We had previously obtained the 
magnetic conditions from a survey of uhe launch area and from measurements at the 
jjosition of the maser in the launcher while the launcher was being erected. The fields 
during liftoff and into the trajectory were obtained from earth field models, and the 
effect of the spinup of the probe at the appropriate eaitli field was incoiporated in the 
simulation. The simulation of the rapid ascent into the vacuum of space was also 
included at the appropriate time in the sequence. 
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The magnetic field encountered during the flight has been determined from standard 
earth models. We found that, only the axial component appears to be important. The 
algorithm developed for the magnetic influence on the probe output frequency has two 
branches owing to hysteresis effects in the magnetic shields. The following algorithm 
was developed in the magnetic calibrations; 

^ = -5.8992 X 10~^^ ®axial , 

~ = -12.7213 X 10"^^ ®axial descending . (19) 

The sign convention is that the plus direction of field enters the earth at its north 
magnetic pole. 

The effect of temperature variations within the maser causing cavity frequency 
shifts was measured via telemetry in terms of the aft-oven heater voltage. The 
frequency variation due to this effect was determined during tests. The following 
relationship for frequency versus aft-oven heater voltage was found: 

^ = -3.60 X 10"^'^ AV . (20) 

The effect of pressure variations measured in the dome enclosing the forward 
assembly of the maser containing the cavity resonator was determined by calibration 
to be as follows: 

Af -12 

-^=-2.26X10 AP (decreasing pressure) , 

-3. 11 X 10~^^ AP (increasing pressure) , (21) 

where AP is in psi. The pressure was measured by telemetry. 

The effect of spin-rate variations on the oulput frequency of the maser was measured; 

^ = (0. 1906 X lO"^*^) - (7.554 X 10“^^) , (22) 

where is in rpm. An optical device on the payload provided rotation-rate data through 
the telemetry system. 


12. EXPERIMENT OPERATIONS 


The final payload weight was 225.49 lb. On May 3, 1976, SAO received an updated 
trajectory prediction from the Vought Corporation for a payload weight of 226. 1 lb. 
Table 1 illustrates the as-flown performance witli respect to the 226. 1-lb trajectory 
prediction of May 3 and with the preflight predictions made from an updated weight 
estimate of 225.49 lb. 


Table 1. Flight performance. 


Parameter 

May 3, 1976 
prediction 
(226. 1 lb) 

Preflight 
prediction 
(225.49 lb) 

Flight 
results 
(225.49 lb) 

Impact of spacecraft 
Latitude, geodetic 
Longitude 

24:5909 N 
5b: 5860 N 

27:5544 N 
48:7257 W 

29:3027 N 
47:4568 N 

Time after liftoff (sec) 

7090.7 

7042.83 

7002.98 

Experiment time (sec) 
(time to x° tracking 
elevation minus 366 
sec) 

7! 6, 6584 sec 
0!7, 6634 sec 


7:9, 6480 sec 
0:4, 6540 sec 

Apogee geocentric radius 

1.6752 X 10^ km 

1.6684 X 10^ 

1.66242 X 10^ km 

Apogee time after liftoff 
(sec) 

3610.87 

3586.84 

3567.21 

Redshift at apogee 

4.306 X 10“^® 

4.295 X 10"^^^ 

4.285 X 10“^® 


The trajectory, apogee altitude, and ejqjeriment time were remarkably close to 
the preflight predictions, and the vehicle performance was completely satisfactory. 


The data-acquisition system at Merritt Island recorded the beat frequency in terms 
of sine and cosine voltages at 14-bit accuracy every 0. 01 sec. The beat data were 
also acquired for "quick-look" purposes by a strip-chart recorder, some features of 
which are reproduced in Figures 18 through 23. 


Figure 18 shows the acquisition of the payload signal while the payload was still 
attached to the fourth stage during the third-stage powered phase of the flight. The 
probe was not rotating, and the signals were being communicated through three para- 
sitic antennas located on the skin of the vehicle and connected to the on-axis anteima 
of the payload by a coupler, which was later disconnected when the payload separated 
from the fourth snage. After third-stage burnout, we can see in the figure the effect 
of the spinup of the payload and fourth stage (while still attached to the third stage) 
before fourth-stage ignition. The phase of the beat signal was broken vq) because of 
the motion of the three separate antennas as the vehicle rotated. 

Figure 19 depicts the end of this broken pattern when the payload separated from 
the bumed-out fourth stage, exposing the dipole antenna. By having located the antenna 
on the spin axis and having made the probe’s exterior surface rotationally symmetrical, 
very little variation of phase with rotation was observed. This figure represents the 
beginning of the useful data. Note that the beat freqi ency was slowing down (time goes 
to the right, with two major divisions equal to 1 sec). 

About 5 min later, the beat signal went through zero, as shown in Figure 20. Here 
we see that the beat has gone from a negative frequency to zero and has reversed in 
sign. At this time, the aspect angle of the on-axis dipole probe antenna as seen from 
the ground was very nearly zero, close to the null of the antenna. The slight misalign- 
ment of the probe's rotation axis with respect to the dipole electrical axis, which 
occurred at payload separation, caused a small but observable modulation of the signal 
amplitude as the probe rotated. Incidental amplitude-to-phase conversion, present in 
all electronic systems, produced the observed phase modulation. Both the rotation rate 
and a modulation of the nutation rate are visible in the trace. These rates agree well 
with those predicted from our measurements both of the moments of inertia and of the 
impulsive torques available from the spinup thrusters. The spin rate also agrees with 
data from an optical spin sensor. 

Figure 21 shows the beats as they continued near apogee about an hour later, and 
in Figure 22, the beat signal is seen going to zero and reversing in sign during payload 
descent. No antenna modulation is visible in Figure 22; nutational motion, as expected, 
had damped out by this time. 
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Figure 19. Sine cosine chart of relativistic effects: 4th stage separation 



Figure 20. Sine cosine chart of relativistic effects: Zero heat on ascent. 





Figure 21. Sine cosine chart of relati\istic effects: Beats near apogee 


h 



P'iKure 22. Sine cosine chart of relativistic effects: Zero beat on descent 




P'ipurc 23. Sine cosine chart of relativistic effects: End of experiment 



The frequency continued to increase as the probe iq>proached the earth, and we 
see the end of the experiment data in Figure 23. The apparent loss of phase coherence 
at about 12 sec before complete loss of signal is probably due to multipath propagation 
at low elevation angles. According to the existing trajectory, the phase apparently 
suffered no loss of continuity down to an elevation angle of 0.4*. 

A time line of the significant events during the mission is shown in Table 2. 

During the mission, we experienced only one loss of phase coherence, which 
resulted from a circuit-breaker dropout in the power supply for the uplink transmitter 
at the Merritt Island station. While this problem was quickly solved and only 1^8^ 
el{q)sed before the probe signals were reacquired, the phase stability of the system 
suffered severely and the data acquired after the dropout are still not understood 
well enough to be included in the results of this e;q)eriment. This problem was further 
complicated by the earlier than anticipated depletion of the liquid ammonia used to 
control the temperature of transponder-translator electronics. After the flight careful 
simulations were made on a similar set of electronics following the same time-tempera- 
ture profiles and only very small phase variations were observed. 

A rough approximation of the shape of the beat frequency versus time is given in 
Figure 24, which is plotted from quick-look frequency measurements. 
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Table 2. Significant mission events, 


Event 

Time (QMT) 

Une3q>ected 

events 


Liftoff 

Firit-atai* burnout 
Flnt-atage aeparation 

11:41:00.24 



Firat-ata(e heat ahield 
ejected 




Second-atas* ignition 
Second-atage burnout 
Second-atage aeparation 




Signal acquired at MILA 
and locked in 

11:43:09.06 



Third-atage Ignition 
Third-atage burnout 
Third-atage apinup atart 
Third-atage qiinup end 
Third-atage aq>aration 
Fourth-atage ignititm 
Fourth-atage burnout 

11:44 



Payload aeparation 

1 1 :46 



Zero beat in redahift 
data on aacending 
branch 

Bandwidth changed from 1 . S to 
0. 5 mHz 

11:49 

12:03 

Payload 

rotation 

rate 

increaaed 


Apogee 

12:40:27 




13:18 

NH 3 ran out 



13:22:47.63 

Loaa of lock; 
MILA failure 



13:23:55.93 

Signal 

reacquired 



13:24:50.86 

MILA exciter 
relocked to 
l^und maaer; 
all lock loops 
closed 

Payload 

rotation 

tale 

constant 

Zero beat in redahift data on 
deacending branch 

13:31 



Firat phaae diacontlnulty; 
liaiA elevation angle 0 . 4' 

13:36:00 
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13. TRACKING DATA AND TRAJECTORY SOLUTIONS 


Tracking data were obtained in the form of d(^pler cycle covints at Bermuda (BDA), 
NTTF (Goddard Space Flight Center. Greenbelt, Maryland), Ascension Island, and 
Merritt Island (MILA), Florida. All stations obtained data in the form of measurements 
of cycles of phase resulting from the change in signal path from Merritt Island, to the 
probe and back to the station. In the case of Merritt Island, which received its own 
uplink transmissions from the probe transponder, Li doppler data represent twice the 
range rate of the probe and are represented as the cycles shown at mixer Ml of 
Figure 8. These data were encoded in the format shown in Figure 16 every 0.2 sec 
on the odd numbered tenths of the second and are interleaved between the data taken 
from mixer M2, which are recorded at even tenths. The other stations recorded doppler 
and telemetry data in the same format, however they did not take one-way data from 
the probe; the three-way data were encoded every 0. 1 sec. 

Our plan was to use NASA's normal data reduction procedure to process the 
doppler data from all stations by fitting them to a ballistic model of the trajectory. 

This procedure evaluates the fit and adjusts the initial conditions of position and 
velocity (state vector) to obtain an c^timum fit to the ballistic model. Corrections for 
tropospheric and ionospheric effects, light- time delay (retardation), residual drag 
from the earth's troposphere, and light pressure are incluueu in the ballistic model. 
During the process of reducing the telemetry data to convert them to useful measure- 
ments, we observed that the rotation rate of the probe changed significantly during 
the mission and that this angular acceleration stopped at a time coinciding with the 
depletion of the ammonia boiler. Careful research of the question of how the ammonia 
was vented from the probe led us to believe that the vent nozzle may have slightly 
misdirected the plume of spent ammonia and that a small residual thrust existed that 
caused angular acceleration of the probe. The existence of a component of thrust 
was present along the probe's spin axis caused a great deal of difficulty in the solu- 
tion of the probe's trajectory from the doppler data. The fact that the uplink tracking 
signal was interrupted v/as also a hindrance to obtaining a clear-cut solution using 
standard techniques. 
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With the usual 20-20 hindsight we now see clearly that a few stretches of data 
from the standard USB range measurement system would have provided extremely 
valuable bench marks of position to verify the solutions obtained by fitting doppler 
data to a dynamical model. 

The first attempts at trajectory solutions by the Goddard Space Flight Center 
showed that some nongravitational forces were at work on the payload and gave some 
insights to tlie possible problems to be encountered in modelling these forces to 
obtain an accurate trajectory using the dynamical methods. A close study of the sus- 
pected thrust from the ammonia boiler observed in the angular acceleration of the pay- 
load gave a time profile of the tangential component. Assuming that the nozzle con- 
figuration remained constant during the ammonia boil off we can infer from this time 
profile what the time profile will be for the axial component. The axial component 
will alter the course of the normally ballistic trajectory the radial component is aver- 
aged out owing to the 115 rpm (approx.) rotation rate. 

This complication led to a reassessment of the problem and further dynamical 
solution methods were attempted at GSFC. Concurrently with these efforts, a non- 
dynamical method was developed at SAG, which amoimted to a strictly geometrical 
intersection of three spheres whose radii were measured in terms of the doppler cycle 
wavelengths counted at the MILA, BDA, and NTTF stations. 

This method of solution is not original and has been long used for tracking vehicles 
under powered flight. Its application depends on establishing a point of departure, or 
some known position during the flight at a known time. This is used to determine 
the number of doppler cycle counts at each of the three stations simultaneously with 
the time of the known position. Here, as mentioned earlier, the use of the USB range 
determining capability would have been very valuable. Fortunately we have data 
from a ‘‘ourth tracking station at Ascension Island, which we used as a constraint on 
the choice of starting coordinates for the geometrical solution. 

By observing the time behavior of the predicted minus the observed residual doppler 
cycle counts from the Ascension station as a function of starting position variations in 
each of three directions, we look for the starting position that gives a zero time slope to 
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this residual. Except for the vQ)link interruption, the cycle counts are continuous with 
very little Jitter. The Ascension residuals give a very distinct and different time signa- 
ture when the originating point is varied 500 meters in each direction. These three signa- 
tures, when expressed in terms of a cubic equation, supply us with a series of poly- 
nomial coefficients for each Ax, Ay, Az from a standard starting point. We can then 
determine the optimum values of Ax, Ay, and Az that make the residual signature 
closest to a straight line with zero slqjc. 

The observations from Ascension were very limited however, and the probe 
never exceeded 6* elevation angle, which led to a problem of effects of tropospheric 
and ionospheric refraction and the question arose regarding how much of the data could 
be considered reliable in steering the starting point. Certainly near apogee, when there 
was little change of elevation angle with time, these effects would be minimum, 
however the extent to which Ascension data before and after apogee could be used was 
an q)en question. Test solutions were made using varying intervals of data centered 
at ^ogee and the behavior and the stationarity of the coordinates of the starting point 
was determined as a function of interval. Using this method, the estimated uncertainty 
in the coordinates put the starting point within a box approximately 40 meters on a 
side. 


A third solution of the trajectory was made by Dr. J. D. Anderson and Mr. A. Liu 
of JPL using the dynamic approach with an axial force parameter modelled on the 
rotational information. This solution emphasized the fitting of the data taken after the 
ammonia was depleted where no thrust was expected and solved for the magnitude 
of the axial thrust coefficient. The standard deviation of the orbital parameters when 
transformed into speed and radial distance by simple conic equations gives a maximum 
radial error of 12 meters occurring near the end of the flight and a velocity error of 
0.04 cm/sec. The magnitude of the thrust was found to be about 163 d>mes and did 
not deviate in direction more than 5° initial velocity vector. The final analysis showed 

the S-Band phase residuals to be less than 0.8 cycles and the maximum doppler error 

-3 

was 0.5 X 10 Hz. A more complete account of this solution is given in Appendix I. 
The ground track of the trajectory is shown in Figure 25. 


75 





This carefully developed solution from JPL finally resolved the nagging questions 
that had long delayed our ability to reach acceptable conclusions about the experiment. 
Table 3 shows for the JPL, SAO and GSFC solutions the magnitudes of the position 
and velocity differences from mean values established at times 1215, 1240 and 1322 
GMT. 


We conclude that the combined error of the position and velocity uncertainties 

-15 

across the intctyal 1215 to 1322 is on the order of ±3.33 X 10 when expressed in 
terms of Af/f given in equation (15). 

The predicted phase accumulation for f = 2117.7 Hz using equation 15 is shown 
in Figure 26. 
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Table 3. 

Time 

1215 


1240 

(apogee) 


Effect of differences between JPL, SAO, and GSFC trajectory solutions on 
redshift prediction. Comparison are shown at 1215, 1240 (apogee) and 1322 
in terms of differences from average altitude and velocity. 

Average Radius from Earth Center, r 
(r = 14,972,895.88 m) 

6 (o 1 from average 
\c^ / 

’'JPL -y- ® 

igAo “ ^ “ +137. 63 m 

^GSFC-^=- 

RMS 6^ 

+1.58 X 10“^® 
-2.7 XIO’^® 
+1.16 X 10"^® 

= 1.93X10"^® 

Average Relative Velocity, v 
(v = 228, 510. 0 cm/sec) 

6(-^ j from average 
\2c/ 

Vjpj^ -v= 0.1 cm/sec 
VgAo - V = -0.7 cm/sec 
VQgFc cWsec 

RMS 6^- 

+0.2 XIO"^® 
-1.8 XlO"^® 
+1.3 XIO"^® 

= 1.29X10“^® 

2c / 

Average Radius from Earth Center 
r= 16,624,474.35 m 

6 from average 

- r = - 77.68 m 
'‘SAO +134. 17 m 

*^GSFC ■ “ 56.48 m 

j 

+ 1.25 X 10"^^ 
-2. 16 X 10"^^ 
+0.9 XIO"^^ 

RMS 6(- 

= 1.53X10"^^ 

V = 70, 710.3 cm/ sec 

6(-^ ) from average 
\2c‘^/ 

Vjpj^ - V = - 2.33 cm/ sec 
VgAo ~ cm/ sec 

'^GSFC"'^~~ 2.33cm/sec 

RMS b(- 

-1.8 Xio"^® 
+3.6 Xio"^® 
-1.8 XIO"^*^ 

= 2.55 X 10"*® 


2c 
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Table 3. (Continued) 


Time 

1322 


r= 11,930,118.6 m 


""JPL 72.16m 

2 

"■sag 122.02 m 

-3 


_1 

RMS = 2 

\c^ / 

V = 449, 864 . 6 cm/ sec 

<i) 

Vjpj^ - v = - 0.7 cm/sec 

-3 

VsAo - V = + 1.2 cm/sec 

6 

Vqs PC - V = - 0. 5 cm/ sec 

l2 

A . 


RMS 


«)■ 


Taking RSS of End Points 1215 to 1322 RMS Values: 



\c^ / 

6 ) 

1215 

±1.93 X 10"^® 

±1.29 X 10“^® 

1322 

±2.68 XiO’^® 

±4.16 X 10"^® 

Overall 

±3.3 XlO'^® 

±4.36 X 10"^® 


Error from Trajectory: ±3.33 X 10 


from average 

.2 XlO"^® 

.8 XlO"^® 

.5 X10~^^ 

.68 X 10“^® 

from average 

.5 XIO'^® 

.0 XIO"^® 

.0 X10“^^ 

.16 X 10"^® 
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UNITS OF 10 RADIANS 



Figure 26. Predicted phase accumulation. 
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14. CALTORATION DATA 


During the 2-hr time span from launch to impact the probe encountered a variety 
of environmental changes. The signal levels in the three microwave links varied widely. 
The effects of these changes on the payload were calibrated during a series of thermal 
vacuum tests. In the caie of magnetic and barometric variations, hystereris effects 
were observed and, where applicable, the correction algorithms have two branches 
for use if the direction of the stimulus is changed. 

The effects of zero g on the payload maser were estimated by upending it and 

-12 

measuring its output frequency shift. This was found to be about 2. 57 X 10 and the 

corresponding cavity shift was found to be 119.3 Hz. One half of this shift, 59.65 Hz, 

was taken as the shift to zero g, and, before launch the cavity frequency was offset 

by this amount so that when it was in a weightless condition it would be close to the 

properly tuned frequency. While our confidence in the structural integrity of the cavity 

was good as it had survived some very severe shock and vibration tests, we developed 

further corrections to be used in the event that the cavity shifted several hundred Hz 

combined with the occurrence of a substantial change in atomic hydrogen flux as 

observed by a variation in maser output power. The cavity shift was measured from 

the redshift prediction using a 500-sec segment of data near apogee. The difference 

-13 

in the output frequency of the maser from the predicted value was -2.91X10 in 

Af/t. From this we deduce that the cavity was low in frequency by 20 Hz. This small 

shift combined with the negligible change in maser ou^ut power level rendered the 

- 13 

correction for flux variation completely negligible. The value - 2.91 X lO in Af/f 
was taken as a fixed frequency offset throughout the mission and was incorporated as 
a constant rate of phase accumulation in the phase prediction as shown in the data 
processing diagram given as Figure 17. 

The only effect that was not simulated during the test sequence was the sustained 
thrust acceleration of up to 20 g from the rocket. Mechanical stress analyses and 
stress testing of components gave us confidence that the probe would survive and 
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that whatever momentary deformations occurred, such as the foam supports of mag- 
netic shields in the maser, would return to their original state within a few minutes. 

A number of effects were identified that produced significant frequency shifts 
and frequency correction formulas were developed during the test phases of the pay- 
load. These are shown in Table 4 "Sources of Possible Bias Errors. " In particular 
we must ke^ in mind the possible errors involved in making these corrections since 
they will contribute uncertainty in the measuring accuracy of the overall experimental 
apparatus. The table shows the magnitude of the variation various perturbations in 
the interval 1215/1322 GMT and the predicted possible bias error estimated at the lo- 
level of confidence. 

The time signatures of these effects is given in Appendix II under two categories: 
Maser Frequency Variation Corrections and System Phase Variation Corrections. The 
latter are very small compared to the frequency effects and their contributions to 
overall bias errors are not considered to be significant. 

The combination of all the corrections for phase and frequency effects is shown 
in Figure 27 in terms of phase versus time. 

Not included in the above is the calibration of the fixed offset error of the probe 
maser versus the ground maser. Owing to changes in the design of the maser for 

space use its temperature and magnetic fields are different from the ground masers. 

-12 

The net offset is 5. 5 X 10 in Af/(, the probe running at a higher frequency by this 
amount. 
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Table 4. Sources 


Parameter 


Axial magnetic 
field from 
trajectory 


00 Maser oven 

temperature 
(oven heater voltage) 

Maser outer can 
pressure 

Probe rotation rate 


T rajectory 
uncertainties 


Correction formula 

= (-5. 8992 ± 0. 567) 

X lO'^^/Oe 
to apogee 

(-12.7213 ± 0.894) 

X 10" ^^/Oe 
from apogee 

= (-3.60 ± 0.60) X 10"* Vv 

(-310.8 ± 4.3) X 10"*Vpsi 

H1.435± 0.48) 

+ (0. 148 ± 0. 0084)nj 
X 10~*Vrpm 

See Table 3 


Root-sum-squares: Estimate of overall possible bias 


possible bias erroic. 


Range of variable 


1215-1240 

1240-1322 

l<r contribution to Af/f 

-0,98 X 10~- Oe 

+3.36 X 10"^ Oe 

±0.6 X 10"*® 



±3. OX 10"*® 

-0.25 V 

+0.25 V 

±1.5X10"*® 

±0. 02 psi 

+0.02 psi 

±1.7 X 10"*® 

+1.4 rpm 

+1.8 ipm 

±4.6X10"*® 


±3.3 X 10"*® 
±6.8 X 10"*® 



RADIANS 




4 


2 


0 



EXPERIMENT TIME (HR^MIN. GMT) 

Figure 27. Combination of all corrections. 



IS. EXPERIMENT RESULTS 


The relativistic effects we are testing, have been discussed in Section 7 and are 
contained in equation (IS): 

♦p(t) - 

The trajectory data and the map of the geopotential provide the information we require 
to generate a predicted value [Af(t)/f]^ as shown in the data-reduction flow plan (Figure 
17). The predictions of the relativistic effects and the systematic phase and frequency 
shifts are all expressed in terms of accumulated phase shifts at the 2117. 7-MHz uplink 
frequency. For convenience these data are all normalized to zero at apogee and have 
been normalized again to zero after the interruption at 1322:47 GMT. Similarly, 
the sine cosine observations given in Figures (18) and (19) from the dc^pler cancelling 
system have been expressed in terms of accumulated phase normalized to zero at 
apogee. 

♦ 

The uncorrected residuals resulting from the difference A^t) = 6 (t) - <t> (t) 

P ® 

are shown in Figure 28. The average slope in the interval 1235:27 to 1245:27 is used 

to determine the probe maser offset frequency owing to launch and zero g effects. 

-13 

The frequency offset, as mentioned earlier, was 3. 1 X 10 on Af/f, which, when 
removed from the plot shown in Figure 28 gives the result shown in Figure 29. The 
predicted systematic phase variations relating to the maser and system shown in 
Figure 27, when removed from the data, give the overall phase residuals shown in 
Figure 30. These have been smoothed by use of a 100-sec filter weighted by a Hanning 
function. The time derivative of the phase, expressed as Af/f, is shown in Figure 3 1 . 
Here we observe most vividly the discontinuity that resulted when the system was 
relocked after the uplink transmitter dropped out at 1322:47. Despite repeated efforts 
to understand this shift, it remains a mystery. All phase lock loops were relocked 

★ _12 
This prediction also includes the 5. 5 X 10 offset between probe and ground masers 

measured before the flight. 
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Figure 28. Redshift phase residuals— uncorrected. 
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Figure 29. Redshift phase residuals- 3. 1 X 10 slope removed. 
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correctly and all monitored quantities were found to be normal. The opportunities for 
conjecture are almost boundless. Because of this uncertainty we have decided not to 
use the data after 1322:47. 

The stabilization of the payload maser frequency after the thrust acceleration and 

mechanical shocks during launch was not simulated in ground tests and we could only 

estimate the time required. From the phase residual plot (Figure 30), we see that this 

restabilization required about 25 min, which was longer than anticipated. At 1215. 

-11 2 

the transient has died out and while there is only 2. 94 X lO change in A<j>/c iq) to 
apogee at 1240:27 these data are retained in the statistical consideration of the results. 
From apogee to 1322 we have a change of 1. 13 X 10~^^ and the e^erimental data set is 
chiefly over this range. 


A statistical description of the data is shown in the regression analysis of fre- 
quency measurements made over separate 100-sec intervals over the time interval. 
Figure 32 shows the least squares fit of a straight line to the residual frequency data. 

The formal value of the quantity e describing the agreement of our data with the redshift 

—6 

hypothesis is e = (5 ± 126) X lO” and we have 


Af 

T\ 


= [1 +(5± 126)X10"®]^ . 

formal c 


We must combine the value of € with the estimate of uncertainty from possible bias 
in the overall system shown in Table 4: 


M 

f 


bias 



, 6.8X10~^^ 

* -10 
1.13X10 


= ±60.2X10“® . 


Combining the bias uncertainty with the formal statistical uncertainty as root-sum- 
squares we have 

^ = (1 ± (5 ± 140) X 10"®] ^ 

c 


as our final result. 
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16. CONCLUSIONS 

The agreement of the experimental data with redshift prediction discussed in the 
previous section depends on two major a priori assumptions. First, we have assiuned 
that the velocity of light is constant and isotropic for distances over which there is a 
substantial difference in gravitational acceleration and potential. Secondly, we have 
assumed that the time dilation, or second-order doppler effect, has been tested e^q>eri- 
mentally at sufficient accuracy to support the redshift conclusion and, further that the 
effect, which to date has been tested only under local conditions, can safely be applied 
to the nonlocal situation where there is substantial curvature in the spacetime metric. 

The question of the precision of the verification of time dilation has been 

* 

reviewed by Newman et al . (1978). To date the smallest ejqjerimental vmcertainty 
has been obtained in the measurement of muon lifetimes by Bailey et al . (1977), who 

o _i/2 

cite lifetimes at different values of Yj “ (1 - Pj ) and conclude that 

^ t ~7 v" \ ^ X 10 or 0. 1% precision. 

if we attribute all the uncertainty in our experiment to the second-order doppler effect 
and assume that the redshift is exactly predicted we have 


— = 1 + (5 ± 140) X 10‘® 
A(v^/2c^) 


and a corresponding precision of 0.014%. This precision is the same as that quoted 

2 2 

for the redshift because the range of v /2c is almost exactly the same as the range 
2 

of aVc since the probe was in nearly free fall. 


The question of a possible asymmetry in the one-way and two-way velocity of light 
measurement can also be addressed if we refer to the relativistically corrected 
* 

We are grateful to Prof. N. F. Ramsey for pointing out that the best test of special 
relativity to date is at the 10~3 level and that our previous conclusions should be 
reconsidered. 


! 
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doppler frequency shift e^qjressions given in equations (11) and (12), and write °12 
as the velocity of light from earth to probe and Cg^ for the velocity in the opposite 
sense. Under these conditions, if we expand equations (11) and (12) we have the 
following: 


f' -f. 



'<>2 - ♦a 




C 


i" - t. 





The output from mixer M3 in Figure 3 gives 


Af 

( 

IT 


i.ff. 

0 0 \ 


0 ■ 2 


■ Af ^ ^23 ^ ^12 ^23 

^ °23 ^®12 ^°23 


^ Af r /_1 1_ 

^ V12 °23/ 


Here we have written r^g = - rgg = r . Writing Ac = Cgg - c^g we have 


Afc i. ^ , r Ac 
^0 


f f. 2 c 


If we ascribe all the observed departure of frequency to the possible anisotropy of the 


velocity of light we conclude that the value of Ac/c is on the order of 3 X lO' 
data set from 1215 to 1322 GMT. 


in the 


The concept of separating the several different effects is not really legitimate since 
the theory we are testing will fail if any of these effects is not correctly observed. What 
is actually being tested is the self consistency of the General Theory of Relativity. 

* 

A recent result on the local invariance of the length of a laser cavity has been reported 
by Brillet and Hall (1979). 
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The measi - e of this self consistency is the e:q)erimentally determined one-sigma 
probability limits of the departure of the frequency shift from the prediction divided 
by the magnitude of the total frequency shift. From this point of view we have 

11 + (2.5 ± 63) X 10'®] 1"^+ 


— 15 

as formal error to which we add the possible apparatus bias of ±6.8 X 10 in Af/f, 

— fi 

which causes a 29.4 X lo" probable error in the determination of 



The conclusion we reach for the self-consistency of the hypothesis is 


•y= [1 + (2.5 ± 70) X 10'®] 



The ejqseriment confirms the self-consistency of the General Theory of Relativity 
at the 70 parts per million level. To date the best test of the Special Theory has been 
made with high-energy muons from CERN particle accelerator and the confirmation 
reported by Bailey et al . (1977) is at the 1000 parts per million level. The previous 
test of the gravitational redshJft by R.V. Pound et al . (1965) using the MSssbauer 
radiation from Fe^^ over a height of 75 feet is at the 10, 000 parts per million level. 
We believe that the experiment described in this report has substantially improved 
the confidence we have in the General Theory of Relativity and the meshing of the 
Special Theory and the Principle of Equivalence upon which it is based. 
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ENCXNEERING MEMORANDUM 
314-187 
7 M*y 1979 


TO; 

FROM: 

SUBJECT: 


T. H. Thorncoa 

A. S. Liu/J. 0. And«rson/D. L. Cain/P. A. Laing 

Gravitational Radahlft Proba Orbit Datamlnatlon from 
Counted Cyclaa Utlog SAO Hydrogan Maaara 


REFERENCE: Gravitational Radahlft Spaca-Proba Exparlmant, R.F.C. Vassot, 

M. W. Levina, GP-A Projact Raport from Smithsonian Astro- 
physical Observatory, April 1977 


I. INTRODUCTION 

% 

An experiment was undertaken In 1976 by Dr. R.F.C. Vassot and 

associates of the Smithsonian Aatrophysical Observatory (SAO) to test 

Einstein's theory of the gravitational redshift effect (Ref. 1). This 

experl:uent used ground based as well ae tocket-bome hydrogen maser 

atomic clock. By comparing thp transmitted probe frequency (controlled 

by the probe hydrogen maser) , with the ground reference frequency it 

was possible to measure the Einstein redshift effect to an accuracy of 
14 

1 part (10 ), which was 100 times more precise than ever achieved. To 
assist SAO In improving this experiment, we produced and sent to SAO a 
trajectory ephemerlu This trajectory was determined from simultaneous 
doppler cracking from 4 stations, 3 of which were equipped with S.AO 
masers. The standard deviations cf the orbital parameters were derived 
from a weighted least square fit to Che doppler data. These standard 
deviations, were trans'^orracd Into speed and radial distance errors by 
simple conic equations. The maximum radial distance error was 12 meters, 
occurring near the end of the flight, and the velocity error was .04 cm/.sec. 


. a 
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With this trajectory, It will ba possible to separate the redshlft effect 
from doppler motion to the accuracy of 1 parc/10^^. 

II. ANALYSIS 

The Cravitatlonal Redshlft Space-Probe was launched In 1976 by a 
Scout D rocket from Wallops Station. The rocket, which reached an altitude 
of 10,000 km, was tracked In a 2-way mode from Metric Island Launch Area 
(MILA), and simultaneously in a 3-way mode by Bermuda, National Test and 
Training Facility (NTTF) , and Ascension Island. 

All stations, with Che exception t;f Ascension Island, were equipped 
with SAO hyrdogen masers. A-priori station locations, given to us by 
Dr. VessoC, and post-priori, station locations determined by our orbit 
determination program, POEAS, are listed in Table I. Coordinates are 
given in a cylindrical system, (distance km from spin axis, CU, distance 
km from equatorial plane, CV, and longitude in degrees, LO). The post- 
priori statistics for all stations were about 7 meters- For the most 
part, station locations were adjusted from two to five meters, however, 
vith the singular exception for longitudes of Bermuda and NTTF. Bermuua's 
longitude .roved eastward by 27 meters and NTTF moved westward by 31 meters 
corresponding to about 1" error in station locations. 
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POST-PUORI VS 

A-PRIORI TRACKING STATION 

LOCATIONS 


CU (KM) 

CV (KM) 

LO (DEG) 

MIU* 

5609.0450398 

3026.0982923 

.379.30656250 


5609.0448033 

3026.1008317 

279.30651987 

BERMUL..^ 

5393.2967222 

3393.4044091 . 

295.34189244 

BERMUDA- 

5393.3020999 

3993.4021880 

295.34213250 

NTTF* 

4963.4455187 

3992.1940424 

283.15744095 

UTTF-^ 

4963.4450874 

3992.2093687 

283.15715810 

ASCEN.IS* 

6517.7157105 

-876.90508990 

345.67282925 

asce:<ms- 

6317.7158143 

-876.90509810 

345.67282940 


* ■ Post- priori value 
“ “ A-priori value 

The dynamical motion of the probe Included percurbacional forces 
caused by a 15th degree and order terrestrial gravitational harmonics, 
point source effects of the moon, sun and all the planets, plus a small 
thrust due to venting from the ammonia boiler. We were able to account 
for this force by first observing the probe's spin ra*te (approximately 
115 rpm) , which changed in time during Che major portion of Che flight. 
Using this spin information we constructed a cubic polynomial repre- 
senc..ng. the time behavior of the spin rata, w. 
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1) w(radians/sec) - .0288 + 5.47 x lO'^t + 1.329 x lO"^ -1.725 x lO'^c^ 
where c is the time Interval In minutes. 

From this, we derived an expression to account for the probe \ gas leak 

acceleration as: 

2) f(km/sec^) - f /ra (1+4.8572641 x 10”^t - 9.4589428 x lO'^t^) 

0 

f ■ 0 after 13:18:00 

where is a constant (in d3mes) to be determined from the tracking data. 

is normalized for a probe mass of 100 kilograms. The ammonia gas was 
exhausted after 13:13:00 because no further change In the spin rate was 
observed. 

In addition to the constant of the acceleration profile, two other 

parameters, U^, U^, describing the direction cosines of the thrust vector, 

were also included and were determined from the data. A third component 

U of this unit vector U is formed by; 

2 

3) U » (1 - U ^ - U 

z X y 

The components of U are expressed in the Earth's mean Equator and 

Equinox of 1950.0 coordinate system. Table II lists the resulting values 

from a combined least square solution of the orbital elements and 

station positions and the thrust vector. The thrust magnitude, f , was 

0 

found to be about 163 d\>nes acting on a 100 kilogram probe. The thrust 
direction did not deviate by more than 5® from the initial velocity 
vector because the "tip-off" direction during separation is controlled 
to about 5°. Using the doppler data, we found the angular separation 
between the velocity vector and the thrust vector to be 5® + 

® C '1- 
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Indicating Chat insofar as the data would allow, the thrust vector and 
the initial velocity vector were aligned at separation. Since the probe 
was spinning at 115 rpm, the spin direction was assumed to be fixed 
throughout this experiment. 

Phase coherent counter reading of the 2-way, 3-way data at each 
second was given us by SAO. These data were averaged to 1 minute doppler 
samples and processed by our double precision orbit program POEAS. Our 
final iteration by POEAS used the counter readings directly and a doppler 
phase or a differential "range" data trajectory was cbnstructed. 

The final analysis showed the S-band phase residuals to be less than 
.3 cycles. The maximum doppler error was 5 mlllihertz, (excluding last 
3 points which were affected by the troposhpere) , or a maximum velocity 
error of .3 ram/sec. A power loss of the uplink transmitter at MILA for 
about two minutes at 13:22:47 caused a momentary loss of data from ail 
cracking stations. Data from Bermuda from 12:07 to 12:15 were also 
rejected. The data quality otherwise was superb. Doppler residuals 
from all stations were on the order of a few mlllihertz. Plots of the 
P0E.\S phase residuals are shown in Figures 1 and 2. 
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f TABLE II 

I THRUST FROM THE 

AMMONIA BOILER VENT 


- .589635 + .066375 

U - .469563 ^ .021204 

y - 

Force on 100 KG Probe 163.177 ^ 1.366 (Dynes) 


Data residuals from Che ocher three staclons, Cracking In a 3-way mode, are 

shown in Figures 3, 4, S. Plots are expressed in mm/sec vs cine, where 1 mm/ sec • 

.015 hz. Bermuda residuals are all on Che order of 5 millihercz (Fig. 3). N'TTF 

e.xhibiced noisier data Chan Bermuda for reasons unknown at present (Fig. 4). 

.Ascension Island's data were received at low elevation angles (below 5”) and at 

Che beginning and end of crack showed large tropospheric effects (Fig. 3). 
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DOPPLER GANCE LUNG SYSTEM PERFORMANCE STATISTICS 

The effeotiveneas of the dqn>ier cancelling system can be illustrated by observing 
the short-term fluctuations in the one-way and two-way doppler data s^arately as they 
are generated in mixers Ml and M2 of Figure 8. By confuting the residual phase 
signature from the predicted cycle counts for the one-way and two-way data, using the 
trajectory we can remove the gross first-order doppler effects and look at the r^id 
fluctuations. Figure m-1 shows the one-way and two-way residuals plotted versus 
time. These data residuals include predictions of the relativistic effects and lig^time 
refraction but do not include the known systematio effects given in Figure 27. By 
digitally reconstructing the function of the analog doppler cancelling system and 
including the corrections we obtain very nearly the same signature in the residual 
phase as shown in Figure 30. It is worth noting that, in this case, the one- and two- 
way data are not obtained at the same time. As mentioned earlier the recording of 
accumulated phase in the one- and two-way channels alternated every 0.2 sec. During 
this interval the spinning probe's orientation chan^d one third of a turn. Evidence of a 
slight asymmetry of the antenna pattern about the spin axis is shown in Figures 18 and 19 
where the null combination of redshift and second-order d(^pler shows the rotational 
and nutational effects. In the nearly simultaneous cancellation from the analog system, 
where the cancellation was limited only by the light time between the earth and probe 
(r < 0.03 sec), this was almost completely removed. However for the 0. 1-sec delay 
between adjacent one- and two-way readii^s, the rotation effect is not cancelled, and 
is the main reason for the sawtooth behavior in the phase signature. 

Some Insight may be available in the residuals after the uplink interruption. We 
notice that both the one- and two-way data show rather violent behavior, which is well 
cancelled before the break but appears to be in^erfectly cancelled after. Propagation 
effects are not likely to be the cause. However, wherever the signal paths are func- 
tionally separated in the transmitter and receiver system there is an opportunity for 
incomplete cancellation since only the transmitter was interrupted and thermal effects 
in the antenna feed beyond the phase sampling point could occur. Further investigation 
in this might prove useful. 






Statistical r^resentati'>ns of the uncanoelled one-a^ data, the cancelled data, 
and the data obtaiTOd between the two reference oscillators at the earth station are of 
Interest, partly to evaluate the effectiveness of the cancellation method, and to assure 
us that the frequmcy stability of the ground-based oscillators was as expected This 
information is contained in the Allan variance plots shown in Figure m-2, which were 
made using data taken during the interval 1215-1322 GMT for the probe versus the 
ground station control maser, P7 (VLG-11 S/N7) and for P7 vs. P6 (VLG-11 S/N6) 
directly at the ground station. 

The statistics for both comparisons are made in the same way; the number of 
samples used is shown with the statistic. Relative frequency measurements between 
the ground masers were made at 0. 83-sec intervals while the probe vs. ground maser 
data are taken in terms of phase measurements every second. Comparisons for time 
intervals less than 0. 83 sec are not possible, however ample data have previously 
been taken between ground masers for shorter time intervals, as shown in Figure 5, 
which describes the performance of ground masers. In Figure 5 the statistics are 
given for measiirements made between exactly similar maser oscillators and their 
performance, given by the variance, is attributed equally to each maser by dividing 
the comparison variance by V%. The noise bandwidth of the measuring system is 1 Hz 
for the probe comparison whereas it is 6 Hz for the ground masers. This has little 
effect for averaging times longer than 100 sec, however for averaging times of about 
I sec the additive white phase noise in the maser receiver system dominates and the 
contribution of this noise to the variance varies as the square root of the noise band- 
width. The behavior expected at 1 Hz-noise bandwidth for the P6-P7 comparison is 
shown dotted. This line is a factor v2/6 below the level shown in Figure 5 at t = I sec. 

From Figure III-2 we conclude that the doppler cancelling system contributed very 
little noise and that the maser instability is the dominant source of inaccuracy in the 
experiment. 

A study of the noise spectral characteristics in the one-way data is given in the 
following memorandum from J.D. Anderson and A.S. Liu of the Jet Prcq)ulsion 
Laboratories. 
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2. Probe maser vs. P7 (ground maser) and P6 vs. P7 (ground maser c<Hnpari8aa) during same time 
interval. Allan variance data taken during fl^it from 1215 to 1322 GMT, June 18, 1976. 
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315.5.212 
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TO: Distribution 

FROM: J. D. A&d«rson/A. S. Liu 

SUBJECT: Noise Characteristics of Oneway S>Band Data from the Gravitational 

Redshlft Space-Probe Experlaent 


Stnanary 

One-way and two-way Dcppler data froa the 1976 gravitational redshlft space-probe 
expsrlnent have been sent to us by R. F. C. Vessot of the Soithsonlan Astro- 
physical Observatory (SAO) . Our detenslnatlon of the trajectory of the Scout D 
rocket trajectory fron the two-way phase data have been described elsewhere 
(Liu et al, EM 314-187, 7 May 1979). This asno describes sooe "quick- look" 
analysis of the one-way data. Power spectra In phase and frequency are 
generated by fast Fourier transform (FFT), and the Allan variance for the 
phase residuals Is cosqiuted. For Foturler frequencies less than lO'^ Hz we find 
that the power spectrum Increases with lower frequency, but because of the 
limited length of the data record (64 min) It Is not possible to determine the 
precise powerlaw dependence on frequency. If a Kolmogorov spectrum Is assumed, the 
phase spectrum obeys the following law. 




An Independent calculation of the Allan variance from one-way phase residuals yields 




5.3 X 10 

“ 7 “ 


-12 


-2 

For Fourier frequencies greater than 1*^ Hz the spectrum Is representative of 
white phase noise (S. ■ 0.07 radians^/Hz) . The spectrum of fractional frequency 
variations can be fit quite well by the following law. 


Sy(f) - (6.1 X 10"^^) f^ (f > 10“^ Hz) 
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and tha aaaoelatad Allan varlanea ia glvan by 


Oy(T) 


1.5 X 10 


-11 


This la About ona ordar of aatultuda too larga for hydrogan aasar ayttan nolaa, 
and nay ba tha raault of tha aliasing of tha payload spin (% 116 rpn) in tha 
apactrua, or it could ba causad by high fraquaney anvironaantal nolaa in tha 
tropoaphara or in tha payload. Wa plan to invastigata this problaa in tha 
future. 

Tha Kolaogorov apactrua ia roughly of tha aaaa aiaa aa that for tha Vllcing S*band 
data at opposition (Araatrong, Eatabrook, and Wahlquist, Sac. 326 lOM, April 
17, 1979). Hovavar, bacausa Viking has an S/X band radio systaa, it is posslbla 
to attribute aoat of tha low frequency nolaa to plaaaa turbulence. In tha case 
of the radahift spaca-prbba, only S*^and data are avallabla, and so it la lapOB<> 
slbla to separata troposphere and ionosphara affects. Tha apactrua ia aada up 
of aoaa uidcnown coahination of ionospheric and tropos^aric turbulence. Hovavar, 
an upper botmd of Oy(1000 a) ■ 1.7 x 1(T^^ can ba placed on tropospheric nolaa 
in the ona>vay space probe radio link, which has a hydrogan aasar on both ends. 
Conditions at Goldstona, where tha air is dryer and where elevation anglaa arc 
higher (radahift data arc on tha ordar of 59*’), ahould result in considerably 
lass tropospheric noise than indicated by our upper bound frms the space probe. 


Power Spactrua 

First, tha bast-fit trajectory fr<» tlM two-way phase data was passed through the 
one-way phase data recorded at a 60 s saapla. Than a quadratic polynoaial was 
fit to the resulting one-way rasiduala to reaova trends causad by a frequency 
bias and a frequency drift in the data. The resulting residuals are shown in 
Figure 1. Next, the first 64 residuals ware subsdttad to an FFT and tha power 
spectrum was calculated. In addition 128 one-way phase data at a one second 
sample near the beginning of the trajectory ware analysed in the same way. The 
coahlnad hl^ fraquaney (1/128 < f < %) and low frequency (1/3840 < f < 1/120) 
spcctnm from the two sets of residuals is slMwn in Flgtire 2. It is our intention 
to compute a spectrum from a long record of ona second data in the future. 

Tha spectrum can be fit reasonably wall by tha following modal 


S^(f) - 263 + 0.070 (radians^/Hs) 


The constant term in this modal represents white phase noise. The standard 
deviation for tha phase noise can ba evaluated from the formula 
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J %<« " 


and for tha uhita phaaa eo«ponant» 


0^ • (.070/2)** ■ 0.19 radians 


Tha 128 ona>sacond phaaa rasiduala vara also augaantad by ona point and then 
128 fractional fraquancy rasiduala vara coaputad by tha foraula 
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* 27T(2.2 X lOO 


Tha powar spactrua was coaputad for thasa rasiduala as wall (?lgura 3) , and it 
was found that tha resulting S^(f) could ba raprasantad by tha following law. 


Sy(f) - (6.1 X 10"^^)f^ 


Tha standard deviation a 


y 


can ba found froa 


and 




1/3(6. I I 
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0 ■ 1.6 X 10 
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Than is a aaall inco^tlbility batwaan o. ** eoaputad froa the two 
apactra, baeauaa for whita phaaa noiaa, with no eorralatlon batwaan sao^las, 
tha followint ralation ahould hold. 



2 0 . 


(2n v^) 


2 


If wa ccwputa froa 


■ hf . wa find that 


S 


♦ 


(2. V,)* 

I 



and for ■ %. • 2.2 x 10®, h ■ 6.1 x 10*^^, 

HO 


S. • 0.0A9 
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which is nasonably close to tha high fnquancy spactrw of Figun 2. Tna 
incoapatibility is not significant. 


Allan Variance 


The Allan variance is a two-sanpla varianca for avaraga fractional fraquancy 
variations ovar a cycla count tias t. Whan we are daallng with phase residuals 
with trends raaovad, and tha count tiaa ia equal to tha saapla tiaa T, tha 
avaraga fractional fraquancy variation ovar T is: 


y,j(T) 


*k+l ~ 

2n V T 
o 


fdxich is si^)ly Doppler fnquancy nsiduals divided by v . 

o 

Wa can also cosqnita y|^ for T > T, but to avoid cornlatlon between saaples of 
they snist be exactly adjacent in tiaa. This can be accoapllshad by defining 
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yj^(2T) - h + y2it.i<T)] 


This produ.:es a Hat of yi^(2T) which la at wat one half the size of the Hat of 
yu(T). Similarly we can produce a list of yj^(4T) which la at beat h the size of 
the list of , and so on until there are too few values of yj^Cl) to compute 

a meanlngfxil estimate of the two-aa&f)le variance. 

The two-sample variance is a little strange, but it Is used In practice^to 
characterize oscillator stability. Normally from a sample of N values y^, we 
would compute the r by 



Then the unbiased estimate of the variance would be given by 


N 


-v2 


^ E (y,-?) 

k"l 


If we have only two sauries of then the above two equations can be combined 
to yield 



The Allan variance Is defined as the average of the estimates of the two-sample 
variance. Therefore, for a list of m values of yj^, the Allan variance Is given 
by: 


a ‘(T) - 

y 


m-1 


0-1 

E 

k-1 




We have applied this formula to the 93 phase residuals at a one minute sample 
rate (T - 60s) . First we computed the differenced Doppler residual and divided 
by 2.2 X 10^ to produce a list of 92 values of yj^CT). Then 46 values of y^^ (2T) were 
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produced end so on. The resulting Allen veriences ere given in the following 
table. 


T 

m 

0 (t) 



y 

60 s 

92 

7.21 X 10“ 
3.87 X 10“,, 
3.38 X 10“f, 
2.75 X 10“,, 

1.80 X 10“'^'’ 

120 s 

46 

2A0 s 
480 8 
960 8 

23 

11 
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The five values of can be represented by the following law. 


a (t) 

y 


5.30 X 10 



,-12 


This is an alternative description of noise to the more familiar power spectrum, 
but in fact they are related through the relation 


IN 


(,) . 2 / s (f) it 

^ V 0 m K ^ 


(irfr) 


For white phase noise (S^ * hf ) this expression can be Integrated to yield 

9 3hf . 

(t) - ^ + r [sin(4TTf t) -•8sin(2irf„T)] 

y (2 ttt)^ 2(2,tt)^ L N N J 


and 


(3hf„) 


T > >1 


For the white phase spectrum of the one-second data, the formula yields 


o 

y 


,_v "V 1.5 X 10 

(X) 
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-14 

The Allan variance for hydrogen Baser systems typically reaches 10 for 
integration tisws of 100 s. Therefore the phase noise in the one-second 
data that we have analysed is an order of magnitude too large for system 
noise. It must be caused by soamthing else. 
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